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ABSTRACT 


This  thesis  explores  the  recent  rise  of  skiing  and  snowboarding  injuries  from 
jumps  fabricated  with  no  scientific  design  process.  It  summarizes  a  previous  method  to 
develop  a  first  order  ordinary  differential  equation  (ODE)  for  the  landing  surface  shape, 
based  on  kinematics  and  dynamics,  which  limits  the  equivalent  fall  height  (EFH)  on 
landing.  These  are  compared  with  theoretical  expressions  for  EFH  in  tabletop  jumps, 
which  are  shown  to  have  linearly  increasing  and  possibly  large  EFH  values  near  the  ends 
of  the  tabletop  and  linear  landing  surface  portions.  Finding  solutions  to  the  ODE  is 
explored,  with  a  large  emphasis  on  determining  the  singular  point  where  the  ODE 
numerical  integration  can  begin.  Analysis  is  conducted  to  determine  a  good  way  to 
design  a  curved  in-run  transition  portion  of  the  jump  that  limits  the  additional  centripetal 
acceleration  on  a  particle  undergoing  the  required  velocity  change  of  direction  of  a  given 
amount.  This  turn  can  be  accomplished  using  a  unique  curve  known  as  a  clothoid,  which 
minimizes  jerk  along  the  path.  The  final  topics  include  a  plan  for  providing  maintenance 
for  safe  terrain  park  ski  jumps,  and  an  algorithm  that  will  assist  a  manager  in  planning 
where  and  how  to  build  a  jump. 
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EXECUTIVE  SUMMARY 


This  thesis  explores  the  rise  of  skiing  and  snowboarding  injuries  from  jumps 
fabricated  with  no  scientific  design  process,  and  presents  a  method  to  design  a  much  safer 
jump.  Recent  skiing  progression  has  resulted  in  a  much  more  aggressive  and  acrobatic 
style  of  skiing.  Current  tabletop  terrain  park  jumps  are  not  designed  with  any  scientific 
process  that  mitigates  injuries.  Furthermore,  a  skier’s  equivalent  fall  height  (EFH) 
increases  as  the  jump  size  increases.  Clearly,  there  is  a  need  for  a  safe  jump  transition, 
ramp,  and  landing  surface. 

A  safe  jump  landing  surface  can  be  designed  through  exploring  the  kinematic 
requirements  for  a  safe  landing  scenario,  where  the  flight  path  angle  is  equal  to  the 
landing  surface  angle.  On  a  tabletop  jump,  angle  equality  like  this  rarely  occurs.  By 
choosing  a  reasonable  and  low  EFH,  it  is  possible  to  limit  the  perpendicular  velocity 
upon  landing.  As  a  result,  a  low  impulse  on  landing  occurs,  which  results  in  a  softer 
landing.  Exploring  these  requirements  leads  to  a  first  order  ordinary  differential  equation 
(ODE)  for  the  landing  surface  shape.  The  EEH  in  tabletop  jumps  have  linearly  increasing 
and  possibly  large  EEH  values  near  the  ends  of  the  tabletop  and  linear  landing  slope 
portions.  Conversely,  the  EEH  in  a  safe  jump  design  remains  constant  at  the  desired 
level. 

The  ODE  has  an  infinite  family  of  solutions  for  given  values  of  the  EEH,  takeoff 
angle  Oq  and  initial  velocity  Vg-  A  skier  flight  path  intersects  these  ODE  solutions  at 
points  that  have  values  of  EEH  of  a  desired  quantity.  The  ODE  has  a  “singular  point,”  or 
starting  point,  where  the  ODE  numerical  integration  can  begin.  This  unique  singular 
point  is  shown  to  lie  along  a  rotated  ellipse,  which  is  generated  for  a  set  of  parameters 
00  and  EEH.  Without  knowing  the  singular  point,  the  only  method  to  integrate  the  ODE 
is  through  estimating  a  conservative  integration  starting  point.  With  the  singular  point 
location,  numerical  integration  becomes  much  more  exact.  Physically,  the  singular  point 
exists  where  the  skier  flight  path  is  perpendicular  to  the  landing  surface. 


Further  analysis  is  conducted  to  determine  a  good  way  to  design  a  curved  in-run 
transition  portion  of  the  jump  that  limits  the  additional  centripetal  acceleration  on  a  skier 
undergoing  the  required  velocity  change  of  direction  of  a  given  amount.  A  circular 
transition,  found  on  design  for  some  Olympic  ski  jumps,  does  not  limit  the  jerk  that  a 
skier  feels  upon  entering  and  exiting  the  transition.  A  better  shape,  which  limits  jerk,  can 
be  accomplished  using  a  unique  shape  known  as  a  clothoid.  A  clothoid  has  a  variable 
radius  along  its  length,  and  it  provides  the  smoothest  link  between  a  straight  line  and  a 
circular  curve.  On  a  clothoid,  the  distance  from  the  start  of  a  transition  and  the  radius  at 
that  point  are  inversely  proportional  to  each  other.  For  a  given  set  of  incoming  and 
outgoing  angles  in  a  transition  region,  an  exact  clothoid  shape  can  be  found.  Comparing 
the  centripetal  accelerations  found  in  a  circular  transition  versus  a  clothoidal  transition 
shows  that  the  clothoid  shape  is  a  better  design. 

Through  careful  sight  designation  and  thoughtful  analysis,  a  ski  area  manager  can 
design  a  safe  in-run,  transition,  takeoff  ramp,  and  landing  surface.  Furthermore,  safety 
and  skier  enjoyment  can  be  continued  by  ensuring  that  jump  maintenance  is  conducted 
daily. 
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I.  BACKGROUND  ON  SKI  JUMPING 


A.  AN  ABBREVIATED  HISTORY  OF  SKIING 

Skiing  is  one  of  the  oldest  sports  in  human  history.  The  earliest  known  accounts 
of  skiing  occurred  over  4000  years  ago,  when  people  in  Scandinavia  skied  for  travelling 
and  hunting  [1].  Proof  of  early  skiing  as  a  method  of  travel  and  hunting  is  documented  in 
Figure  1,  which  depicts  a  petroglyph  of  a  skier  discovered  in  Norway  [2].  More  recently, 
in  the  past  150  years,  skiing  has  transitioned  from  a  utilitarian  sport  into  one  of  pleasure, 
adventure,  and  recreation.  As  skiing  gained  popularity,  ski  clubs  and  competitions  began 
appearing  in  Europe  around  1850,  and  later  in  the  United  States  [3].  Skiing  culture 
continued  to  grow  as  equipment  improved,  allowing  skiers  to  perform  better  and  enjoy 
themselves  more  than  ever  before  [4] . 


Figure  1.  Earliest  known  ski  drawing,  from  [2] 


B.  THE  RISE  OF  FREESTYLE  SKIING 

The  broader  sport  of  skiing  has  many  different  aspects  that  have  evolved 
throughout  the  past  few  hundred  years.  Some  of  these  smaller  niche  styles  include  mogul 
skiing,  racing,  ski  ballet,  jibbing,  powder  skiing,  aerials,  ski  jumping,  skiercross,  and 
others. 
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1.  Olympic  Jumps  and  Aerials 

Ski  jumping  gained  popularity  in  Norway  in  the  mid  1800s,  and  has  beeome  a 
facet  of  the  Winter  Olympics  and  many  other  competitions.  In  this  aspect  of  skiing, 
skiers  slide  down  a  path  and  glide  off  a  jump.  They  fly  through  the  air  for  as  long  as 
possible  before  landing.  Skiers  are  rated  for  their  style  and  how  far  they  go  before 
landing  [5].  Figure  2  shows  a  ski  jump  in  Einsiedeln,  Switzerland. 


Figure  2.  Ski  jump  in  Einsiedeln,  Switzerland,  from  [5] 


Aerial  ski  jumping  is  another  kind  of  skiing.  In  this  sport,  skiers  launch  from  a 
nearly  vertical  jump  and  perform  many  flips  and  twists  before  landing  on  a  steep  landing 
slope.  Judges  grade  each  jump  on  the  complexity  of  the  aerial  maneuvers  that  skiers 
perform  in  the  air,  and  how  good  the  landing  is  [6] .  Eigure  3  shows  a  skier  performing  an 
aerial  trick  off  a  jump. 
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Figure  3.  A  skier  performing  an  aerial  trick,  from  [7] 


2.  Steep  Skiing 

Since  the  mid-1990s,  skiing  has  progressed  in  a  highly  acrobatic  and  aggressive 
direction.  Along  with  the  growth  in  popularity  of  the  action  sports  industry,  such  as 
motor  bike  racing,  skateboarding,  surfing,  wakeboarding  and  other  sports,  skiers  and 
snowboarders  quickly  learned  to  ride  on  steeper  terrain  and  launch  off  bigger  jumps  than 
ever  before.  Rather  than  being  confined  to  a  set  of  rules  on  how  to  ski,  the  “freeskiing” 
movement  continues  to  grow  every  year.  Today,  skiers  are  pushing  the  boundaries  of 
where  it  is  possible  to  ski.  Skiing  in  locations  and  conditions  that  were  once  considered 
unthinkable  and  impossible  just  a  few  years  ago  is  now  quite  attainable.  In  Figure  4,  a 
skier  (in  the  circle)  is  contemplating  his  next  turn  down  an  extremely  steep  mountain  in 
Alaska  [8]. 
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Figure  4.  Steep  skiing  in  Alaska,  from  [8] 


3.  Terrain  Parks 

Within  the  past  15  years,  man-made  terrain  parks  grew  from  simple  jumps  to  a 
series  of  obstacles,  including  features  such  as  large  jumps,  halfpipes,  quarterpipes,  and 
rails  to  slide  on  [9] .  This  unprecedented  change  in  a  ski  area’s  terrain  led  to  a  new  way  to 
think  of  what  is  possible  while  riding  on  skis  or  a  snowboard. 

a.  Components  of  a  Current  Terrain  Park  Jump 

The  current  industry  standard  for  a  terrain  park  jump  consists  of  shaping 
snow  into  a  “tabletop”  jump  [10].  Figure  5  presents  the  five  main  components  to  a 
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conventional  tabletop  terrain  park  jump:  the  in-run,  the  transition,  the  takeoff  ramp,  the 
tabletop,  and  the  landing  slope  where  skiers  exit  the  jump  area.  There  is  also  a  location 
known  as  the  “knuckle,”  which  is  single  spot  where  the  tabletop  changes  to  the  landing 
surface. 


Figure  5.  Profile  for  a  typical  terrain  park  jump 

b.  Non-existent  “Design  ”  of  a  Table-Top  Jump 

Recently,  two  journals  provided  analysis  of  terrain  park  jumps  and 
methodologies  for  improvement.  Hubbard  [10]  and  McNeil  and  McNeil  [11]  discussed 
the  fabrication  of  today’s  tabletop  jumps  and  provided  guidelines  to  make  them  safer. 
Despite  these  publications,  most  ski  areas  do  not  use  any  engineering  methods  in  the 
design  and  construction  of  terrain  park  ski  jumps.  Currently,  ski  areas  privately  envision 
and  control  the  construction  of  tabletop  jumps.  This  process  varies  by  ski  area,  and  is 
managed  solely  based  on  the  experience  of  the  ski  area  staff.  Most  often,  the  jump  shape 
and  dimensions  are  based  on  past  jumps  that  the  staff  has  seen,  built,  or  imagined.  As  a 
result,  there  is  an  ad-hoc  method  in  building  a  jump  ostensibly  to  limit  injury  while  still 
providing  a  fun  feature  for  skiers  to  jump  off  [10].  Most  often,  after  construction  of  the 
jump  with  a  snow  cat  and  shovels,  a  few  brave  skiers  test  its  suitability  by  skiing  off  it. 
Depending  on  the  landing  and  experience  of  the  skier,  the  jump  builder  might  slightly 
adjust  the  jump  to  make  it  better.  This  process  may  continue,  as  the  skiers  provide 
feedback  to  the  jump  builder,  who  in  turn  adjusts  the  shape  of  the  jump,  including  the 
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takeoff  angle,  length  of  the  flat  tabletop,  and  length  and  slope  of  the  landing  area.  The 
end  result  is  a  tabletop  jump  that  is  then  opened  to  the  public  for  any  skier  to  enjoy,  while 
accepting  the  risks  involved. 

On  a  tabletop  jump,  a  skier  hopefully  approaches  the  jump  at  a  fast  enough 
speed  so  that  he  will  fly  through  the  air  past  the  flat/horizontal  tabletop  portion,  and  land 
on  his  skis  somewhere  along  the  linear  landing  area.  The  experienced  skier  often  knows 
how  fast  to  be  skiing  for  a  given  jump,  but  an  inexperienced  skier  may  find  himself  out  of 
place.  Too  little  speed  at  the  beginning  of  the  jump  may  result  in  a  landing  on  the  flat 
tabletop.  Conversely,  too  much  speed  at  the  beginning  of  the  jump  may  result  in  landing 
deep  into  the  ending  of  the  landing  surface,  which  typically  becomes  flatter  as  the  feature 
ends.  The  severity  of  impact  (which  can  result  in  injuries)  depends  on  the  velocity  of  the 
skier  [10].  These  facts  are  the  risks  that  accompany  current  terrain  park  ski  jumps. 
Figure  6  shows  a  typical  tabletop  ski  jump  found  at  most  major  ski  areas,  and  Figure  7 
shows  a  skier  performing  an  aerial  trick  after  leaving  the  takeoff  ramp.  Notice  the 
horizontal  tabletop  and  the  linear  landing  slope. 

■ - -  J 


Figure  6.  Current  tabletop  ski  jump,  from  [12] 
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Figure  7.  A  skier  jumping  off  of  a  terrain  park  jump,  from  [13] 


c.  Severe  Injuries  on  Table-Top  Jumps 

As  terrain  park  ski  jumps  have  become  more  popular  for  skiers  to  perform 
tricks  off,  serious  spinal  injuries  have  become  more  prevalent.  The  progression  also  led 
to  an  increase  in  skiing  injuries,  especially  to  the  neck  and  head  [14].  Although  there  are 
highly  successful  ski  safety  programs  and  equipment,  the  injuries  continue.  The 
fabrication  of  these  terrain  park  features  themselves  may  have  contributed  to  the  injuries 
sustained  by  the  skiers  and  snowboarders  [11]. 

d.  Scientific  Method  Needed  to  Design  Terrain  Park  Jumps 

The  severity  of  the  injuries  that  some  have  sustained  while  on  terrain  park 
jumps  motivates  a  need  for  a  scientific  method  to  build  terrain  park  jumps.  Although 
there  are  many  terrain  park  features,  such  as  rails,  boxes,  and  halfpipes,  this  thesis  will 
focus  solely  on  man-made  terrain  park  jumps.  It  will  explain  a  procedure  to  design  safer 
ski  jumps  to  minimize  skier  injuries,  and  will  explain  the  mathematical  conditions  that 
lead  to  the  design  of  both  an  in-run  transition  and  a  landing  surface. 
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For  continuity,  consider  a  reference  system  with  the  skier  travelling  from 
left  to  right.  A  skier  will  ski  down  an  initially  linear  surface,  enter  a  curved  transition 
region  and  then  jump  or  glide  off  a  takeoff  ramp.  The  origin  of  the  (x,y)  coordinate 
system  is  at  the  jump  take  off  point,  which  is  the  last  point  of  contact  with  the  takeoff 
ramp’s  snow  surface.  The  x  axis  is  horizontal.  All  measurements  of  distance  will  be  in 
meters,  angles  will  be  in  both  degrees  (in  written  explanations),  and  radians  (in 
calculations).  Additionally,  all  references  in  this  paper  to  a  skier  are  to  men,  such  as 
"he/his/him."  This  is  for  simplicity;  it  does  not  imply  that  all  skiers  are  men.  Anywhere 
that  "he/his/him"  is  mentioned,  it  can  certainly  be  replaced  with  "she/hers/her".  Also,  all 
references  to  "ski/skier/skiing"  can  certainly  be  replaced  with 
"snowboard/snowboarder/snowboarding".  This  analysis  applies  for  all  jumping  winter 
sports  athletes,  men  and  women,  skiers  and  snowboarders. 

In  keeping  with  previous  research,  and  for  simplicity  in  modeling  the 
flight,  this  thesis  will  ignore  rotations  and  tricks  of  the  skier  while  in  the  air,  and  follow 
the  motion  of  the  skier  center  of  mass  [10].  Additionally,  the  term  “safe”  is  meant  to  be 
used  in  the  context  of  a  skier  landing  with  his  skis  against  the  ground.  This  study 
assumes  that  a  skier  will  be  able  to  land  properly;  it  is  not  a  study  concerning  a  faulty 
landing. 
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II.  DESIGN  METHOD  EOR  SAFE  SKI  JUMP  LANDING 

SURFACES 


A.  SAFER  JUMP  LANDINGS  THROUGH  IMPULSE  MINIMIZATION 

As  a  skier  leaves  the  takeoff  ramp  and  flies  through  the  air  at  a  certain  velocity 
and  takeoff  angle,  there  will  be  a  unique  point  that  he  will  impact  the  landing  surface. 
After  leaving  the  takeoff  ramp,  but  prior  to  impact,  when  the  skier  is  in  the  air,  the 
velocity  of  the  skier  is  most  likely  not  parallel  to  the  landing  slope  [10].  After  landing 
impact,  all  of  the  velocity  is  in  the  tangent  plane  of  the  landing  surface.  As  a  result,  the 
perpendicular  velocity  Vj_  instantaneously  changes  from  nonzero  magnitude  for  flight 
velocity  v^,  to  zero,  with  sliding  velocity  Vj,. 


Figure  8.  Perpendicular  velocity  changes  upon  landing  impact 

The  concept  of  force  is  understood  by  Newton’s  Second  Law  as  equivalent  to  the 
product  of  mass  and  acceleration  [15].  Similarly,  impulse  is  the  integral  of  force  with 
respect  to  time  [16].  Upon  landing,  the  skier  will  feel  a  force,  and  therefore  an  impulse, 
of  the  landing  surface  against  his  skis. 

Force  =  F  =  ma  (2.1) 

Impulse  ~  ^  ~  j  ~  j  ~  '^'^b  ~  '^'^a  (2-2) 
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The  skier  velocity  before  impact  v^,  has  a  perpendicular  component,  while  the  velocity 
after  impact  Vi,  does  not.  As  a  result,  the  impulse  is  proportional  to  v_i  by  a  mass 
constant  m. 

/  =  rnvj_  (2.3) 

In  order  to  make  the  impulse  small,  the  perpendicular  velocity  must  be  small.  By 
limiting  Vj_  on  a  jump  landing,  the  impulse  is  limited,  and  the  skier  will  experience  a 
safer  landing  than  if  the  impulse  is  large.  Even  if  a  skier  flexes  his  legs  to  absorb  the 
shock  of  landing,  he  will  still  be  subjected  to  the  same  impulse. 

A  skier  will  have  a  lower  chance  of  sustaining  an  injury  if  his  velocity  upon 
landing  is  nearly  parallel  to  the  landing  surface.  The  more  nearly  parallel  to  the  landing 
slope  that  a  skier’s  landing  velocity  is,  the  less  perpendicular  the  velocity  will  be  against 
the  landing  surface.  This  ideal  situation  results  in  a  small  normal  impulse,  and  the  skier’s 
chance  of  injury  will  be  lower  [10].  As  the  landing  slope  angle  cp  and  skier  trajectory 
angle  9  become  increasingly  similar  in  magnitude,  the  skier  will  land  more  safely.  As  a 
result,  the  chance  of  injury  will  be  less,  especially  if  the  skier  fails  to  land  on  his  skis. 

Because  limiting  is  not  instinctively  understood,  an  analogous  quantity  can  be 
used,  known  as  equivalent  fall  height  (EFH).  Equivalent  fall  height  is  the  height  that 
results  in  velocity  v^in  a  Ig  environment.  The  measure  of  EFH  is  known  through  the 
energy  equation  h  =  I2g.  The  essential  requirement  is  that  is  made  to  correspond 
to  a  reasonable  EFH  [10].  Figure  9  shows  a  skier  flight  path  with  a  safe  landing  surface. 
The  perpendicular  component  of  velocity  is  the  sine  of  the  difference  of  the  safe  landing 
slope  angle  and  the  skier  flight  path  angle. 

Vj_  =  V  sin((p  —  0)  (2.4) 
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When  ^  =  6,  the  perpendicular  component  of  velocity  is  zero.  In  this  ideal 
situation,  there  is  no  perpendicular  velocity;  all  velocity  is  parallel  to  the  landing  surface 
and  no  impulse  is  required  on  landing. 


B.  DETERMINATION  OF  SKIER  VELOCITY  AT  TAKEOFF 


Landing  safely  on  a  surface  is  related  to  the  EFH  and  velocity  upon  landing.  In 
the  safest  landing  condition,  the  flight  path  angle  and  landing  surface  angle  have  the  same 
value  at  impact  (0  =  (p).  This  section  will  explore  the  velocity  of  a  skier  at  the  ramp 
takeoff  point.  Most  explained  concepts  and  equations  in  sections  B  and  C  of  this  chapter 
are  a  summary  of  Hubbard’s  work  [10]. 


The  angle  of  the  flight  path  measured  from  horizontal  0  is  determined  from  the 
vertical  y  and  horizontal  x  velocity  components  of  the  skier. 


6  =  tan  ^ — 


V-. 


(2.5) 


The  angle  (p  measured  from  horizontal  is  found  by  the  inverse  tangent  of  the  rate  of 
change  of  the  landing  surface  y  distance  to  its  x  distance. 


(p  =  tan  ^ 


dys 

dx 


(2.6) 


A  skier  will  leave  the  jump  takeoff  by  either  jumping  off  the  top  edge  or  by 
gliding  smoothly  from  the  takeoff  point.  If  the  skier  simply  glides  off  the  jump,  the 


11 


initial  takeoff  velocity  Vg  is  parallel  to  the  jump  takeoff  angle  Pq,  relative  to  the  horizon. 
These  horizontal  and  vertical  components  of  initial  velocity  are  given  by: 


o 

O 

II 

o 

(2.7) 

Vyg  =  ViSinpo 

(2.8) 

However,  not  every  skier  will  simply  glide  off  the  jump  into  the  air.  He  may 
decide  to  increase  the  velocity  at  takeoff  by  popping  off  the  top  of  the  jump  takeoff  point. 
Nearly  all  of  his  change  in  velocity  will  be  perpendicular  to  the  velocity  vector  by  an 
amount  Av.  In  modeling,  this  velocity  increase  Av  will  have  an  upper  limit  ^j2gK,  where 
K  is  the  vertical  height  the  skier  can  jump  in  meters.  Typical  recreational  skiers,  wearing 
boots  and  skis,  will  jump  a  maximum  of  approximately  0.25  meters.  The  corresponding 
Av  is  therefore  about  2.21  meters  per  second. 


When  the  skier  glides  off  of  the  jump,  the  initial  velocity  Vg  is  parallel  to  the  ramp 
angle,  Pq.  As  a  result,  the  jump  takeoff  angle  Oq  is  equal  to  /?o.  In  this  scenario,  Vq  =  Vi- 
If  the  skier  jumps  with  additional  velocity  Av,  the  takeoff  angle  changes  by  an  additional 

amount  A6  =  tan~^  — .  The  initial  velocity  then  becomes  Vq  = 


With  the  added  boost  from  the  skier,  the  new  takeoff  angle  Gq  is  Gq  =  Pq  +  AG. 
The  faster  a  skier  is  moving,  the  higher  he  will  be  in  the  air  after  leaving  the  jump. 
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C.  DEVELOPMENT  OF  AN  ORDINARY  DIFFERENTIAL  EQUATION  FOR 

A  SAFE  LANDING  SURFACE 

For  a  safe  landing  surface,  the  primary  requirements  are  to  have  skier  impulse 
from  the  surface  both  limited  to  a  reasonable  value,  and  to  be  reasonable  no  matter  where 
a  skier  lands.  The  velocity  components  along  the  path  as  a  function  of  time  are  known 
from  basic  physics  projectile  equations  (neglecting  aerodynamic  forces). 

Vx  =  Vxo  =  Vo  cos  00  (2.9) 

Vy  =  Vyo  =  Vo  sin  00  -  gt  (2.10) 

Integrating  these  expressions  with  respect  to  time  allows  us  to  find  the  (x,y) 
position  of  the  skier  at  any  desired  time. 

X  =  V;,ot 


y  =  Vyot 


(2.11) 

(2.12) 


Solving  for  time  yields: 


t  = 


X 


(2.13) 


Vo  cos  00 

By  eliminating  time,  the  parabolic  vertical  position  y  of  a  skier  flight  path  can  be 
written  in  terms  of  only  the  variable  x  and  parameters  Vo  and  0o. 

,2 


y  =  X  tan  0o  — 


gx‘ 


(2.14) 


2  (Vocos  0o)^ 

Solving  for  Vo  results  in  an  expression  for  the  initial  velocity  a  skier  needs  in 
order  to  reach  a  certain  point  (x,y). 


Vo 


x^g 


1 2  (x  tan  00  —  y)  (cos  0o)^ 


(2.15) 


The  total  velocity  has  v  and  y  components. 


V  =  IV;c^  +  Vy^ 


(2.16) 
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Substituting  Equations  2.9  and  2.10  in  Equation  2.16  results  in  an  equation  for 
skier  velocity  involving  parameters,  t,  Vq  and  Oq. 

V  =  ^J{VQ  cos  GqY  +  (vo  sin  Oq  -  gty  (2.17) 

Substituting  for  t  yields  the  following  equation  for  the  skier  velocity  at  an 
arbitrary  {x,y)  position  in  flight. 


>1 


(2.18) 


Using  Equation  2.15,  Equation  2.18  can  be  rewritten  in  terms  of  x,  y,  and  6q.  This 
gives  the  expression  for  the  velocity  at  any  {x,y)  coordinate  along  the  flight  path. 


V  = 


x^g 


2{x  tan  6q  —  y)  cos^  6q 


-  ^gy 


(2.19) 


The  rate  of  change  in  vertical  position  is  the  derivative  of  Equation  2.14. 


dy  Vy 

—  =  —  =  tan  6  =  tan  9^ 
dx  Vx 


gx 

(VqCos  9o)^ 


(2.20) 


With  Equation  2.20,  it  is  possible  to  solve  for  the  flight  path  angle  9  in  terms  of 
known  parameters  g,  Vq,  and  9o,  and  variable  v. 

1  gx 

9  =  tan  ^(tan9o  -  z - (2.21) 

Cvocos9oy 

Eurther  simplification  (by  substituting  Equation  2.15  into  2.21)  results  in  the 
flight  path  9  involving  only  terms  of  known  variables  x,  y,  and  parameter  9q. 

2v 

9  =  tan~^( - tan9())  (2.22) 

X 

Building  upon  these  established  formulas,  it  is  possible  to  develop  a  differential 
equation  that  must  be  satisfied  by  a  safe  surface  landing  shape  y^.  The  safest  place  to 
land  is  where  the  landing  slope  and  the  skier  flight  path  angles  are  equal,  and  there  is  no 

perpendicular  velocity.  Erom  Equation  2.4  and  knowing  that  =  ^j2gh,  it  is  possible 
to  solve  for  <p,  the  landing  surface  angle. 
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(2.23) 


(p  =  sin 


+  e 


Using  Equation  2.21,  the  landing  surface  angle  can  be  found  in  terms  of  known 
parameters  g,  h,  and  Oq  . 


(p 


=  sin-i  +  tan-\tan  6^  -  - - ^ 


\  V  J 


(VqCos  6>o) 


r) 


(2.24) 


With  Equation  2.6,  the  rate  of  change  of  the  landing  surface  angle  ^  can  be 


found. 


dys 

tan  (D  =  ——  =  tan 
dx 


^tan  00 


(VoCOS0o)^J  \  ^  JJ 


(2.25) 


Substituting  the  relationship  for  velocity  (Equation  2.17)  in  Equation  2.25  yields  a 
differential  equation  for  the  safe  landing  surface. 


dys 

dx 


/ 


tan 


tan  ^  I  tan  Oq 


\ 


9X  \ 
(vocos  9oy) 


+  sin 


/ 

\\ 

2gh 

(voCosOq)^  +  (voSinOo 

9X  .2 
Vocos9oy  /  / 

(2.26) 


This  ordinary  differential  equation  (ODE)  for  the  safe  landing  surface  is  in  terms 
of  variables  v,  y,  and  parameters  Vq  ,6o,  g,  and  h.  While  g,  h,  and  6o  are  constant  values 
(with  00  varying  as  desired),  Vq  is  known  already  in  terms  of  x  and  y.  With  the  known 
expression  for  Vq  (Equation  2.15),  simplification  results  in  an  expression  for  the  rate  of 
change  of  the  landing  surface  in  terms  of  variables  x,y,  and  parameters  h  and  0o. 

h  ^ 

^  (2.27) 

4(xtan9o  —  y)  (cos  6o)^  ^ 


,  2y 

dys  _  tan(tan  ^( - tan6o)+sin  ^ 

dx  ^ 
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A  safe  landing  surface  solution  is  found  when  the  skier  flight  path  y  crosses  the 
landing  surface  This  occurs  when  y  =  y^. 


dys 

dx 


2ys 

tan{tan  - tan  Oq)  +  sin 


-1 


^  4  (x  tan  9q  —  y^)  (cos 


ys 


(2.28) 


This  first  order  ODE  must  be  satisfied  in  order  to  allow  for  a  safe  landing  slope 
for  any  given  value  of  Oq  and  desired  value  of  h.  It  is  an  ODE  for  the  dependent 
variable  y^  (the  landing  surface  height)  as  a  function  of  x  (the  dependent  variable).  Note 
that  the  gravity  constant  g  no  longer  appears  in  Equation  2.28. 


The  EEH,  also  denoted  as  h  in  the  equations,  is  a  quantifiable  measure  of  the 
acceptable  distance  that  a  person  feels  comfortable  with  falling  to  the  ground.  Eor 
example,  most  adults  might  be  comfortable  with  an  EEH  of  one  meter,  but  anything  over 
that  is  less  acceptable.  With  this  safe  ski  jump  model,  it  is  possible  to  design  a  safe 
landing  surface  for  a  given  EEH  (preferably  at  around  one  meter).  As  a  result,  a  jump 
can  be  built  so  a  skier  will  “fall”  only  the  EEH  value,  unlike  what  happens  today  in 
existing  tabletop  ski  jumps. 


D.  NUMERICAL  SOLUTIONS  OF  SAFE  LANDING  SURFACE  ODE 

The  solution  to  the  ODE  cannot  be  arrived  at  analytically.  Rather,  it  must  be 
solved  numerically  using  a  program  such  as  MATEAB.  Eurthermore,  integration  must  be 
completed  backwards,  from  the  terminal  conditions.  Because  the  ODE  is  not  Eipschitz 
continuous  (to  be  discussed  in  detail  later),  it  is  not  possible  to  integrate  forward  from  the 
origin.  Eigure  1 1  is  a  plot  of  five  solutions,  corresponding  to  y^  (vertical  landing  surface 
limit)  of  Im,  -Im,  -3.5m,  -7m,  and  -8.5m,  and  a  40  meter  horizontal  jump  x^.  These  five 
solutions  are  chosen  from  an  infinite  family  of  solutions. 


16 


15 


10 


-10  - - - 

-15 1  ^  ^  ^  i  ^  i  ^  i 

5  10  15  20  25  30  35  40 

X  (m) 

Figure  11.  Plot  of  five  safe  landing  surfaces  for  a  40m  jump 

The  faster  a  skier  is  travelling  at  takeoff,  the  further  horizontally  he  will  fly.  In 
order  to  limit  the  landing  impulse,  the  value  of  (p  must  increase  as  the  distance  increases. 
As  a  result,  all  solutions  have  monotonically  decreasing  slope  [10]. 

Overlaying  skier  flight  paths  (from  Equation  2.14)  for  different  takeoff  velocities 
onto  safe  landing  surfaces  shows  how  these  surfaces  can  result  in  the  same  EFH  for  all 
velocities.  Figure  12  shows  the  same  safe  landing  surfaces  in  figure  11,  but  also  includes 
skier  flight  paths  from  velocities  of  5  to  40  mph,  in  5  mph  increments.  The  intersection 
of  every  landing  spot  and  the  flight  paths  all  have  the  same  EFH,  no  matter  which  jumper 
path  one  is  on.  In  this  example,  EFH  is  arbitrarily  chosen  to  be  Im  and  the  jump  takeoff 
angle  is  25  degrees. 
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Figure  12.  Plot  of  five  illustrative  safe  landing  surfaces  (heavy  lines)  with  eight 

jumper  paths  (light  lines) 

Superimposing  a  conventional  table-top  jump  onto  one  safe  landing  surface  (from 
an  infinite  family),  it  becomes  apparent  that  safe  ski  jump  design  calls  for  a 
monotonically  decreasing  slope  as  a  function  of  v,  while  a  table-top  jump  maintains  a 
generally  linear  landing  slope.  Figure  13  shows  that  depending  on  the  tabletop 
dimensions,  there  is  at  most  one  or  two  landing  points  that  these  two  different  jump 
landing  surfaces  have  in  common.  Clearly,  safe  landing  surfaces  differ  from  what  is 
currently  being  built  at  ski  areas. 
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In  [10],  Hubbard  outlines  an  algorithm  for  safe  landing  surface  design. 

1 .  Decide  what  EFH  is  desired.  A  possible  value  is  1  meter,  but  this  must  be 
chosen  by  the  designer  based  on  safety  considerations. 

2.  Decide  on  entire  jump  size,  so  it  fits  into  a  given  area  (x^,  y^). 

3.  For  a  fixed  take-off  angle  Pq,  choose  a  safe  surface  that  passes  through 
(xiiyi)-  Overshoot  may  occur,  and  limiting  Vi  may  be  needed. 

4.  For  a  free  take-off  angle  Pq,  choose  a  safe  surface  that  passes  through 
(xL.yL),  so  that  all  maximum  reasonable  Vj’s  can  be  accommodated. 
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E.  TABLE-TOP  JUMPS  DO  NOT  LIMIT  EQUIVALENT  FALL  HEIGHT 


As  mentioned  in  Chapter  I,  table-top  jumps  are  often  associated  with  serious 
skiing  injuries.  This  section  will  show  that  on  a  table  top  jump,  the  equivalent  fall  height 
generally  increases  as  the  horizontal  flight  path  distance  x  increases.  These  results  will 
be  contrasted  with  a  safe  jump  design,  and  the  EFH  values  will  be  compared  for  similar 
size  jumps. 


I.  Equivalent  Fall  Height  as  a  Function  of  Distance 

As  previously  stated,  perpendicular  velocity  Vj_  can  be  better  understood  by 
measuring  the  corresponding  EFH,  which  is  related  to  Vj_  by  EFH  =  I2g.  Equation 
2.29  is  formed  by  using  Equation  2.4  as  a  substitute  for  Vj_. 


EFH  = 


sin^((p  —  9) 

2g 


(2.29) 


This  equation  can  be  rewritten  in  terms  of  x  by  substituting  the  expressions  for  v 
(Equation  2.19)  and  9  (Equation  2.22).  Equation  2.30  expresses  EFH  as  a  function  of 
variables  x,y  along  a  flight  path  and  jump  design  parameters  9^  and  (p. 


EFH  (pc')  = 


4  cos^  6a{x  tan  Bq  —  y) 


y 


sin 


2y 

(p  —  (tan  ( - tanOo)) 

X 


(2.30) 


On  a  tabletop  jump,  a  skier  will  have  two  general  locations  where  he  can  land: 
before  or  after  the  knuckle,  which  has  coordinates  (x^,  —yt).  Figures  14  and  15  show  the 
two  cases.  The  length  of  the  tabletop  is  x^,  the  height  of  the  takeoff  point  above  the 
tabletop  isyt,  and  the  landing  surface  angle  is  The  time  at  the  zenith  (top)  of  the 
skier  flight  path  is  denoted  by  t^. 
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If  a  skier  lands  on  the  tabletop  (Figure  14),  a  simple  expression  for  EFH  can  be 
determined.  At  the  zenith  of  the  trajectory,  the  EFH  has  two  portions,  /ij  and  /i2.  In  this 
case,  /12  is  the  height  of  the  takeoff  point  above  the  tabletop,  and  is  the  height  of  the 
flight  path  zenith  above  the  takeoff  point. 


hi 


(Vyo)^ 

2g 


(2.31) 


The  initial  vertical  velocity  is  related  to  the  initial  horizontal  velocity  by  the 
tangent  of  the  takeoff  angle.  Equation  2.32  is  a  rewritten  version  of  Equation  2.5  with 
initial  conditions  =  v^o,  Vy  =  Vyo  ,  and  0  =  Oq,  which  occurs  when  a  skier  boost  off 
the  jump. 


Vyo  =  v^o  00 


(2.32) 


Substituting  Vy^  from  Equation  2.32  into  2.31  results  in  an  expression  for  /i^. 


hi  = 


tan^  00 

2g 


(2.33) 


After  substituting  v^q  =  Vq  cos  0o  (Equation  2.9)  and  simplifying,  can  be 
expressed  in  terms  of  Vq,  Bq,  and  g. 
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Vq^  sin^  6q 
2^ 


(2.34) 


K  = 


Substituting  the  previously  determined  initial  veloeity  Vo  (Equation  2.15)  into 
Equation  2.34  results  in  an  expression  for  in  terms  of  x,  Oq,  and  y^,  so  that  =  — y. 


/ii 


tan^  6q 

4(x  tan  00  +  Jt) 


(2.35) 


The  final  EEH  is  found  by  adding  the  additional  constant  /i2. 


EFH  =  /ii  + 


h2=yt  + 


x^ tan^  00 
4(x  tan  00  +  yt) 


(2.36) 


Equation  2.36  gives  the  EEH  for  a  skier  landing  on  a  tabletop.  The  next  case 
shows  the  EEH  for  a  skier  landing  on  a  downward  slope  (Eigure  15). 


When  a  skier  lands  on  the  downward  sloping  linear  landing  surface,  the  height  y 
is  given  by: 

y  =  — yt  +  (x  —  Xt)tan(p  (2.37) 

Substituting  for  y  (Equation  2.37)  in  Equation  2.30  gives  an  expression  for  EEH 
when  landing  after  the  knuckle  in  terms  of  x  only,  together  with  desired  parameters  9o,  (p, 
Xt,  and  yt- 
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EFH= 


4cos^  dg(xtandg  —  (y^  +  {x  —  X[)tan  (p) 


-Vt 


—  (x  —  Xi)tan(p 


(p  —  (tan 


.i2(y,  +  (x-x,)tan(p 


—  tan  dg) 


(2.38) 


The  two  general  landing  cases  occur  when  the  skier  lands  on  the  tabletop  before 
the  knuckle  (governed  by  Equation  2.36)  or  on  the  downward  sloping  landing  surface 
after  the  knuckle  (Equation  2.38). 


2.  Analysis  of  Tabletop  Equivalent  Fall  Height  Plots 

Equivalent  fall  height  sensitivity  can  be  determined  by  choosing  a  reasonable 
linear  takeoff  angle  &o  and  individually  varying  the  parameters  tabletop  X(-,  jump 
height  yt,  and  landing  surface  <p.  Each  tabletop  jump  landing  scenario  is  researched 
individually. 


a.  Flat  Landing  Equivalent  Fall  Height  Analysis 

When  a  skier  lands  before  the  knuckle  (x  <  x^),  his  EEH  is  governed  by 
Equation  2.36.  Eigure  16  shows  that  changing  y^  to  realistic  values  (Om,  Im,  2m,  3m) 
while  keeping  6q  constant  at  25°,  the  EEH  values  are  spaced  fairly  evenly  and  increase 
roughly  linearly  as  the  horizontal  length  of  the  jump  increases.  The  value  for  (p  in  this 
case  is  zero,  as  the  skier  lands  on  the  flat  tabletop.  Eor  a  large  jump,  the  EEH  values 
become  quite  dangerous. 
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Figure  16.  EFH  of  a  flat  tabletop  landing  as  a  function  of  jump  distance  x  for  several 

values  of  and  constant  9q  =  25° 


Similarly,  for  several  values  of  the  takeoff  angle  Oq  while  holding 
constant  shows  again  that  the  EFH  nearly  increases  linearly  as  the  horizontal  length  of  the 
jump  increases.  Figure  17  shows  EFH  vs.  x  for  several  0o values  and  a  constant  y^  at  one 
meter.  The  knuckle  position  is  arbitrarily  large,  since  the  landing  will  always  be  on 
the  tabletop.  Again,  the  EFH  grows  to  dangerous  levels  as  the  jump  size  increases. 
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Figure  17.  EFH  of  a  flat  tabletop  landing  as  a  function  of  jump  distance  x  for  several 

values  of  0q  and  constant  =  Im 

In  both  cases  of  individually  varying  the  parameters  y^  and  00^  the  EFH  is  shown  to 
increase  linearly  with  jump  distance. 

b.  Angled  Landing  Surface  Equivalent  Fall  Height  Analysis 

When  a  skier  lands  after  the  knuckle  {x  >  Xf),  his  EFH  is  governed  by 
Equation  2.38.  Plots  of  EEH  as  a  function  of  jump  distance  are  shown  in  Eigure  18. 
Here,  y^  varies  from  0  to  3  meters  while  all  other  parameters  are  held  constant  (x^  = 
10m,  00  =  25  ^  =  —30°)  . 
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horizontal  length  of  jump  x  (m) 


Figure  18.  EFH  of  landing  on  a  constant  downward  slope  as  a  function  of  jump 
distance  .r  for  several  values  of  y^and  constant  parameters  X(-=10m, 

00=25°,  and  (p=  -30° 

Every  graph  in  Eigure  18  starts  at  its  respective  height  above  the  tabletop 
and  increases  to  the  value  of  x^.  At  the  knuckle,  the  EEH  dramatically  decreases  to  a 
comfortable  level.  However,  thereafter,  the  EEH  again  increases  linearly  without  bound. 
Choosing  larger  values  of  x^,  while  keeping  other  parameters  constant,  results  in  a  very 
similar  outcome  of  EEH  behavior.  Eigures  19  (x^  =  20m)  and  20  (x^  =  30m)  show  that 
moving  the  knuckle  position  does  not  change  the  overall  behavior  of  the  EEH  on  a 
tabletop  jump.  In  each  case,  the  EEH  increases  as  the  horizontal  jump  length  increases  to 
the  knuckle  position.  Then,  although  the  EEH  drastically  decreases  at  the  knuckle,  it 
immediately  begins  to  increase  again.  Tabletop  jumps  today  are  built  with  generally 
constant  linear  landing  surfaces.  This  type  of  landing  cannot  protect  a  skier  from  a  large 
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EFH  (m) 


EFH  if  he  lands  far  out  horizontally  from  the  takeoff  point.  Clearly,  the  linear  landing 
surface  is  not  an  ideal  design,  as  EFH  is  never  consistently  limited. 


horizontal  length  of  jump  x  (m) 


Figure  19.  EFH  of  landing  on  a  constant  downward  slope  as  a  function  of  jump 
distance  v  for  several  values  of  yjand  constant  parameters  X(-=20m, 

00=25°,  and  (p=  -30° 
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horizontal  length  of  jump  x  (m) 

Figure  20.  EFH  of  landing  on  a  constant  downward  slope  as  a  function  of  jump 
distance  .r  for  several  values  of  and  constant  parameters  X(-=30m, 

00=25°,  and  (p=  -30° 

Figures  18,  19,  and  20  verify  that  varying  the  knuckle  position  Xf  for 
different  y^  and  constant  0o  and  cp  does  not  affect  the  roughly  linear  dependence  of  the 
EFH  and  horizontal  length  of  the  jump  x.  For  jumps  with  no  step-down  to  the  tabletop 
(yt  —  0)’  the  skier  flight  path  is  symmetric  about  the  zenith  of  the  trajectory.  In  this  case, 
the  landing  angle  will  be  equal  in  magnitude  to  the  takeoff  angle.  If  the  magnitude  of  the 
flight  path  angle  0  is  greater  than  the  magnitude  of  the  landing  surface  angle  (p,  a  skier 
will  be  able  to  land  anywhere  on  the  landing  surface.  The  closer  the  landing  is  to  the 
knuckle  (but  after  it),  the  lower  the  EFH  is.  However,  if  the  magnitude  of  the  flight  path 
angle  0  is  less  than  the  magnitude  of  the  landing  surface  angle  (p,  there  is  a  physical 
region  where  landing  becomes  impossible.  This  region  can  be  imagined  by  travelling 
through  the  landing  surface  and  arriving  at  the  top  of  it;  an  impossibility.  Figures  18,  19, 
and  20  all  show  this  impossibility  foryj-  =  0.  In  Figure  18,  this  region  exists  at 
approximately  10  to  12  meters.  After  about  12  meters,  it  is  possible  for  a  skier  to  land  on 
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the  downward  slope,  and  the  EFH  continues  to  increase.  In  Figure  19,  this  region  exists 
at  approximately  20  to  23  meters,  and  in  Figure  20  it  is  around  30  to  34  meters. 
Increasing  can  make  it  possible  to  land  anywhere  on  the  downward  slope. 


horizontal  length  of  jump  x  (m) 


Figure  21 .  EFH  of  landing  on  a  downward  slope  as  a  function  of  jump  distance  x  for 
several  values  of  (p  with  constant  parameters  X(-=10m,  y,.  =  Im,  and 

00=25° 

Figure  21  shows  that  varying  the  landing  surface  angle  cp  while  keeping 
other  parameters  constant  also  results  in  a  generally  increasing  EFH.  From  the  origin  to 
the  knuckle,  all  values  of  EFH  are  the  same,  as  <p  plays  no  part  in  landing  in  this  region. 
For  a  symmetric  flight  path,  where  =  Om,  the  takeoff  angle  0o  equals  the  landing 
angle  (p.  In  this  case  (  0o  =  (p),  the  EFH  is  zero.  For  cases  where  0o  >  (p,  landing  on 
the  downward  surface  is  possible,  and  the  EFH  is  small,  positive,  and  reasonable.  For 
situations  where  6q  <  (p,  the  EFH  becomes  negative,  which  correlates  to  the  impossible 
region  previously  described.  A  rule  of  thumb  in  terrain  park  jumps  today  states  that  the 
angle  of  the  landing  surface  should  equal  the  takeoff  angle.  This  rule  attempts  to  control 
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the  EFH  that  a  skier  experiences  upon  landing.  However,  as  the  recent  figures  show, 
landing  on  the  tabletop  can  result  in  large  and  unsafe  EFH.  Furthermore,  although  the 
EFH  decreases  when  landing  occurs  beyond  the  knuckle,  it  immediately  begins  to 
increase  and  can  be  unsafe  for  large  x  values. 

Comparing  these  EFH  tabletop  plots  for  the  EFH  from  a  similar  sized  safe 
jump  (Figure  22)  presents  some  astounding  differences.  Tabletop  jumps  have  just  a  few 
instances  where  the  EFH  would  be  considered  comfortable  to  a  skier,  while  a  safe  jump  is 
designed  with  the  low  EFH  requirement  first  and  foremost  in  mind.  Any  landing  point 
along  a  safe  landing  surface  has  the  same  EFH  as  any  other.  Figure  22  shows  that  no 
matter  where  a  skier  lands  on  a  safe  surface,  his  EFH  will  be  the  same  as  any  other. 
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Figure  22.  EFH  of  landing  on  a  safe  landing  surface  as  a  function  of  jump  distance  x 

with  constant  parameters  0o=25°  and  yt=lm 

If  a  jump  has  consistent  and  reasonable  EFH  values,  injuries  can  be 
minimized.  Furthermore,  skiers  will  know  that  the  jumps  they  enjoy  are  built  with  safety 
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in  mind.  Equations  2.36  and  2.38  are  analytic  expressions  for  EFH  in  tabletop  jumps. 
Equation  2.36  models  the  EEH  felt  by  a  skier  landing  on  the  flat  tabletop,  while  Equation 
2.38  models  the  EEH  felt  by  a  skier  landing  on  the  linear  downward  sloping  landing 
surface.  Figures  16  and  17  show  the  numerical  implementations  of  Equation  2.36,  and 
Figures  18-21  show  the  numerical  implementations  of  Equation  2.38.  An  overarching 
and  alarming  consistency  arises  in  all  of  the  plots:  the  EFH  increases  roughly  linearly 
with  horizontal  distance  x  no  matter  what  the  ramp  height  is,  no  matter  where  the 
knuckle  position  is,  and  no  matter  what  the  angled  landing  slope  ^  is.  If  the  takeoff 
angle  Bq  is  equal  to  (p,  and  y^  =  0,  it  is  possible  to  have  small  EFH  just  past  the  knuckle. 
This  is  sometimes  called  the  “sweet  spot”  [17].  However,  for  all  landings  past  the 
knuckle,  the  EFH  increases  roughly  linearly  as  the  distance  jumped  increases.  A  tabletop 
jump  landing  is  safe  only  if  a  skier  lands  in  the  “sweet  spot”  just  past  the  knuckle. 
However,  this  critical  requirement  puts  much  constraint  on  the  skier,  who  will  not  be  able 
to  greatly  control  his  landing  location  after  takeoff.  Eanding  before  or  after  the  sweet 
spot  is  not  ideal,  as  the  EFH  can  be  at  a  dangerous  level. 

Constant  downward  sloping  landing  surfaces  are  therefore  not  safe  if  the 
landing  distance  x  is  significantly  larger  than  the  knuckle  position  x^,  even  if  cp  =  6 q.  As 
a  result  of  this  analysis,  one  can  conclude  that  large  tabletop  jumps  are  not  safe,  despite 
what  adjustments  one  makes  to  the  knuckle  position,  takeoff  ramp  angle,  and  landing 
surface  angle,  if  the  tabletop  is  long,  or  if  the  landing  position  on  the  downward  surface  is 
large. 
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III.  SINGULAR  POINT  OF  THE  SAFE  SLOPE  ODE 


A.  MOTIVATION  FOR  THE  EXISTENCE  OF  A  SINGULAR  POINT 

In  Figure  11,  it  is  clear  that  the  five  safe  surface  solutions  emanate  from  a 
common  point,  and  all  appear  to  have  the  same  slope  at  that  point.  This  “singular  point” 
occurs  only  at  one  location  along  any  safe  landing  surface.  Physically,  the  singular  point 
must  be  below  and  to  the  right  of  the  origin  (take  off  point).  For  a  given  EFH  and  takeoff 
angle  Bq,  there  is  a  unique  jumper  path  that  hits  all  of  the  safe  surfaces  perpendicular  to 
the  landing  surface.  This  single  skier  jump  path,  with  a  specific  initial  velocity  Vq,  passes 
through  the  singular  point.  Furthermore,  there  is  a  different  singular  point  for  every  value 
of  Vq,  as  the  skier’s  flight  path  change  as  Vq  varies. 

Decreasing  the  initial  velocity  that  a  skier  has  upon  takeoff  decreases  the 
horizontal  distance  at  which  he  will  impact  the  landing  surface.  Continuing  to  decrease 
the  velocity  increases  the  difference  between  the  skier  flight  path  angle  6  and  the  landing 
surface  angle  (p.  The  difference  between  these  two  angles  increases  to  a  limit  ofn/2, 
which  occurs  when  the  skier  velocity  is  perpendicular  to  the  landing  slope. 

This  singular  point  is  the  point  where  forward  integration  can  begin  in  a 
MATLAB  automated  integration  program.  Rather  than  integrating  backwards  only, 
knowing  the  (x,y)  coordinates  for  the  singular  point  for  a  given  EFH  and  6^  enable  the 
ODE  (Equations  2.26,  2.27,  and  2.28)  to  be  solved  forward  numerically  from  the  singular 
point.  Prior  to  knowing  the  singular  point  coordinates,  the  best  method  to  integrate  the 
ODE  is  by  making  a  best  estimate  of  where  the  singular  point  is,  and  integrating  forwards 
from  there,  or  backwards  to  there. 


Rewriting  Equation  2.23  results  in  Equation  3.1,  which  is  an  expression  for  the 
conditions  required  at  the  singular  point.  No  velocity  gets  transmitted  to  a  tangential 
slope  component;  rather  all  of  the  velocity  is  perpendicular  to  the  safe  landing  surface. 


(3.1) 
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This  is  true  when  sin  ^  =  1  =  This  is  equivalent  to  stating  that  the 

numerator  is  equal  to  the  denominator  in  the  argument  of  the  inverse  sine  function  in 
Equation  3.1. 

V  =  ^J2gh  (3.2) 

Substituting  the  available  expression  for  initial  velocity  Vq  (Equation  2.15)  into 
that  for  velocity  along  the  flight  path  v  (Equation  2.18)  results  in  a  new  expression  for  v 
in  terms  of  v,  y,  and  Oq  . 

^  _  x^g  +  gCxtandp  -  ly 

2(xtan0o~y) 

At  the  singular  point,  all  of  the  velocity  is  contained  in  the  perpendicular 
component.  As  a  result,  Equations  3.2  and  3.3  can  be  set  equal  to  one  other. 


Equation  3.4  is  a  relationship  that  the  (x,y)  coordinates  of  the  singular  point  must 
satisfy  for  a  given  h  and  dp  ■  Simplifying  results  in  an  expression  involving  x  and  y  (both 
unknown),  dp  (desired  value  for  a  unique  jump  takeoff  angle,  such  as  25  degrees),  and  h. 
The  result  is  Equation  3.5,  which  has  an  infinite  number  of  solutions,  since  there  is  one 
relationship  between  two  variables.  The  next  sections  will  discuss  a  solution  to  finding 
the  unique  singular  point  for  given  h  and  Op . 

x^  =  (4cos^6p)ih  +  y)(x  tanOp  —  y)  (3.5) 

B.  METHOD  TO  DETERMINE  SINGULAR  POINT 

In  order  to  locate  the  singular  point  for  a  given  h  and  Op,  the  (x,y)  solution  must 
satisfy  Equation  3.5.  One  of  the  solutions  is  the  coordinates  of  the  singular  point. 
However,  the  singular  point  cannot  simply  be  any  ordered  pair  that  satisfies  this  equation. 
It  must  also  possess  the  characteristic  of  a  true  singular  point  at  this  location:  that  the 
skier  flight  path  angle  is  perpendicular  to  the  landing  surface  only  at  this  point.  Plotting 
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Equation  3.5  on  the  (x,y)  plane  results  in  an  ellipse.  This  section  will  make  a  conjecture 
about  where  the  singular  point  lies  on  this  ellipse,  and  then  show  that  the  singular  point  of 
the  safe  landing  surfaces  lies  at  the  intersection  of  the  rotated  ellipse  with  its  semi-minor 
axis. 


1.  Characteristics  of  a  Rotated  Ellipse 

The  general  form  for  any  quadratic  curve  (circle,  ellipse,  hyperbola,  and  parabola) 
is:  [18]. 


-I-  2c2xy  +  2cyx.  -\-  2c^y  -I-  Cg  =  0 

After  defining  the  variables  P  and  Q,  [18] 


P  = 


Cl  C2  C4 

^2  ^3 

C4  C5  Cg 


A  0 


Q  = 


>  0 


p 


(3.6) 


(3.7) 

(3.8) 

(3.9) 


This  quadratic  curve  is  an  ellipse,  where  R  =  +  C2. 

Figure  23  shows  a  general  ellipse  with  its  major  and  minor  axes,  with  the  major 
axis  defined  as  the  longest.  The  semi-major  and  semi-minor  axes  are  half  of  each  full 
length  axis,  respectively.  The  axes  are  perpendicular  where  they  intersect  at  the  center  of 
the  ellipse.  Equation  3.10  gives  a  formula  for  the  length  of  the  semi- major  axis  length  Ug 
in  terms  of  the  coefficients  in  Equation  3.6  [18]. 
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Minor  Axis 


Figure  23.  An  ellipse  showing  major  and  minor  axes 


Cle 


2(CiC5^  +  +  C6C2^  -  2C2C4C5  -  CiCgCfi) 

ic2^  -  4ciC3)(7(ci  -c^y  +4C2^  -  (Ci  +  C3)) 


The  semi- minor  axis  length  bg  can  also  be  found  [18]. 


be 


>1 


+  C3C42  -I-  C6C2^  —  2C2C4C5  —  CjCgCg) 

(C2^  -  4ciC3)(-y(c7^^c^)2T4^  -  (Ci  -f  C3)) 


(3.10) 


(3.11) 


The  center  of  the  ellipse  (Xg,  jg)  is  given  by  the  following  expressions  [18]. 

C3C4  -  C2Cs\ 

-  C1C3  J 

/C1C5  -C2C4 

V  C2^  -  C1C3 


(3.12) 

(3.13) 


If  an  ellipse  is  rotated,  the  counterclockwise  angle  of  rotation  from  horizontal  a 
can  be  defined  as  either: 


1 .  The  angle  through  which  the  ellipse  is  rotated  from  the  horizontal,  which 
is  the  same  as  the  angle  between  the  x-axis  and  the  semi-major  axis,  or 

2.  The  angle  between  the  y-axis  and  the  semi- minor  axis. 
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Figure  24.  Ellipse  rotated  at  angle  a,  showing  the  major  axis  and  minor  axes 
This  counterclockwise  angle  a  is  found  through  one  of  the  following  expressions 

[18]. 


a  =  0  for  C2  =  0  and 


n 


a  =  —  for  C2  =  0  and  >  C3 


tan  i(— — ) 

a  = - - - for  €2^0  and  <  C3 

tan~^(  — ) 

a  =  —  H - ^ - for  C2  ^  0  and  >  C3 


(3.14) 

(3.15) 

(3.16) 


(3.17) 


2.  Specific  Rotated  Ellipse  Equation  in  Terms  of  h  and  Bq 

Combining  like  terms  of  Equation  3.5  reveals  an  equation  for  an  ellipse  rotated 
about  the  origin. 
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(3.18) 


—  4  cos^  6q  tan  Ooxy  +  4  cos^  —  4  cos^  h  tan  6q  x 
+  4h  cos^  00  y  =  0 

Comparing  coefficients  in  Equation  3.18  with  the  general  form  in  Equation  3.6, 
the  values  for  to  Cg  are  easily  identified. 

Cl  =  1 

C2  =  —2  cos^  6q  tan  6q  =  —2  sin  9q  cos  6q 
C3  =  4cos^  00 

C4  =  —2  hcos^  00  tan  6q  =  —2h  sin  6q  cos  6q 
C5  =  2h  cos^  00 
C6  =  0 

Notice  that,  for  this  rotated  ellipse,  C2  is  never  zero  for  any  value  of  0o  ^  0. 
Comparing  coefficients  Ci  and  C3  shows  that  for  reasonable  ski  jump  angles  (less  than  60 
degrees),  Ci  <  C3.  At  0o=  60  degrees,  c^  =  C3.  The  mostly  like  scenario  is  for  any 
takeoff  angle  0o  less  than  60  degrees,  which  is  quite  reasonable,  considering  that  a  skier 
must  be  comfortable  with  jumping  at  this  steep  angle.  Therefore,  we  will  not  consider 
the  cases  where  angles  are  greater  than  60  degrees  (cj  >  C3),  and  thus  we  will  only  be 
using  the  third  case  for  a  (Equation  3.16). 

Equation  3.18  presents  a  condition  h  and  Oq  that  the  singular  point  must  satisfy. 
Using  the  known  values  for  the  coefficients  Ci  to  Cg,  it  is  possible  to  rewrite  Equations 
3.10  and  3.11  in  terms  of  h  and  0o  only  as: 

h^^2  cos  0n  (2  sin  0n  —  1) 

a,  = - ,  °  = -  (3.19) 

{cos^  00  —  1)(Y(1  —  8  cos^  00  +  32  cos^  0o)  —  4  cos^  0o  —  1) 

/1V2  C05  00  (2  sin  Or,  —  1) 

be  = - ,  °  = -  (3.20) 

(cos^  00  —  1)(— ^(1  —  8  cos^  00  +  32  cos"^  Oq)  —  4  cos^  Og  —  1) 

Similarly,  the  center  of  the  ellipse,  found  in  Equations  3.12  and  3.13  can  be 
rewritten  in  terms  of  h  and  0o,  as: 
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(3.21) 


Xg  =  h  tan  Oq 

3.  Rotated  Ellipse  Analysis 

The  singular  point  ellipse  changes  as  h  and  6q  vary.  In  order  to  plot  this  ellipse 
Equation,  it  is  necessary  to  choose  an  angle  6q  and  a  value  for  h  (EFH).  Notice  that  the 
origin  (0,0)  is  a  solution  to  Equation  3.18.  Plotting  it  produces  an  ellipse  rotated 
clockwise  at  some  angle.  Consider  the  following  plots  (Figures  25,  26,  27)  for  various 
reasonable  values  of  h  and  6q  . 

y 


Figure  25.  Ellipse  of  Equation  3.18  for  0o=l(*°  and  h=\ 


1  Plots  such  as  this  one  were  made  using  a  free  online  source,  http://www.wolframalpha.com,  a 
computational  knowledge  engine. 
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y 


Figure  26.  Ellipse  of  Equation  3.18  for  6q=25°  and  h=l  m 


Figure  27.  Ellipse  of  Equation  3.18  for  0o=45°  and  h=\  m 
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Plotting  five  ellipses  together  in  Figure  28  (for  a  constant  h  but  different  values 
of  6q)  verifies  that  the  origin  is  common  to  all,  while  the  ellipse  rotation  angle  increases 
as  6q  increases.  All  of  these  ellipses  contain  the  origin.  The  unique  singular  point  for 
each  lies  somewhere  along  the  edge  of  the  ellipse,  and  although  the  origin  is  always  a 
solution  to  the  ellipse  equation,  this  is  not  the  singular  point.  Rather,  the  origin  of  the 
ellipse  plot  is  the  origin  of  the  jump  takeoff. 


Figure  28.  A  family  of  rotated  ellipses  for  h=l  meter,  0o=  20°,  30°,  40° 

The  skier  jump  takeoff  angle  6q  relates  to  the  angle  of  the  rotated  ellipse.  As 
00  increases,  it  is  clear  that  a  increases,  but  at  what  rate  with  respect  to  0o?  Substituting 
the  known  values  of  c^,  C2,  and  C3  in  for  the  definition  of  a  in  Equation  3.16  yields  the 
following  expression  for  a,  in  terms  of  the  takeoff  angle  0o. 
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a  = 


(3.23) 


tan  ^(■ 


—4  cos  Oq  sin  Oq 
1  —  4cos^  9q 


) 


2 


The  value  of  6q  is  the  only  variable  that  will  change  the  value  of  a.  Notice  that 
gravity  is  not  involved.  Although  h  changes  the  size  of  the  ellipse,  it  does  not  affect  the 
angle  a.  Plotting  Equation  3.23  for  a  general  case  a  and  Oq  results  in  Figure  29. 


In  this  graph,  6q  values  are  along  the  x-axis;  a  values  are  along  the  y-axis.  Both 
are  in  radians.  This  arctangent  curve  appears  to  be  nearly  linear  (although  not  exactly), 
with  an  approximate  slope  of  0.7.  This  plot  verifies  that,  for  a  quick  approximation  of  a, 
simply  multiply  Oq  by  0.7.  That  is,  a  =  0.7 Oq. 

4.  Conjecture  for  Singular  Point  Location 

Where  exactly  will  the  singular  point  be  along  the  ellipse?  As  a  skier  jumps  from 
the  left  to  the  right,  the  singular  point  must  be  at  a  location  below  and  to  the  right  of  the 
origin.  From  inspection  due  to  the  physical  location  of  the  jumper  path  and  safe  slope 
intersection,  it  would  appear  near  the  lower  right  portion  of  the  ellipse.  This  final  section 
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on  rotated  ellipses  will  conjecture  that  the  location  of  the  singular  point  of  a  safe  landing 
surface  lies  at  the  intersection  of  a  rotated  ellipse  with  its  semi-minor  axis,  as  seen  in 
Figure  30. 


Figure  30.  Singular  point  location  on  a  rotated  ellipse 


Equation  2.15  gives  the  initial  velocity  required  to  reach  a  given  (x,y)  point. 
Using  the  conjecture  for  the  location  of  the  singular  point,  it  is  possible  to  solve  for  Vq  in 
terms  of  g,0o,  and  the  singular  point  location  (x^p^y^p).  This  new  expression  is 
Equation  3.24. 


Uo 


^sp  9 


2  (Xsp  tan  00  -  Tsp)  (cos  0o)^ 


(3.24) 


Equation  2.18  gives  the  velocity  at  any  point  along  the  skier  flight  path.  Again, 
using  the  conjecture  of  the  location  of  the  singular  point,  it  is  possible  to  solve  for  v  in 
terms  of  g,  0o,  and  the  singular  point  location  (Xsp,ysp)-  This  new  expression  is 
Equation  3.25. 
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V 


(Vq  cos  OqY  +  (vq  sin  Oq 


9^sp 


Vq  cos  00 


(3.25) 


Table  1  compares  impact  velocities  by  varying  0o,  keeping  h  constant,  and 
assuming  the  conjectured  location  of  the  singular  point.  The  first  velocity  (line  4)  is  the 
velocity  at  impact  using  Equation  3.25.  The  second  velocity  (line  5)  is  the  velocity  at  the 


actual  singular  point,  which  is  defined  to  be  ^j2gh. 


00,  in  degrees 

10 

20 

30 

40 

1 

EEH  (meters) 

1 

1 

1 

1 

2 

Xsp  (conjectured) 

0.235 

0.485 

0.768 

1.114 

3 

jsp  (conjectured) 

-0.986 

-0.940 

-0.847 

-0.670 

4 

Velocity  at  impact  (m/s) 

Erom  Equation  3.25 

4.428 

4.428 

4.428 

4.428 

5 

Velocity  at  true  singular  point  (m/s) 

Erom  V  =  ^2gh 

4.428 

4.428 

4.428 

4.428 

Table  1.  Test  results  to  find  velocity  at  singular  point  for  constant  h  and  increasing  6q 


Notice  that  the  two  velocity  rows  are  identical  for  each  value  of0o.  This  is 
because  if  a  skier  lands  at  the  conjectured  singular  point,  it  does  not  matter  what  angle  he 
jumps  from.  His  velocity  at  impact,  found  in  Equation  3.25,  is  the  same  as  the  freefall 
velocity  (line  5).  These  results  seem  to  verify  the  conjecture  that  the  singular  point  of  a 
safe  landing  surface  lies  at  the  intersection  of  an  associated  rotated  ellipse  with  its  semi¬ 
minor  axis.  However,  there  is  one  more  compelling  reason  why  this  location  is  correct. 

One  characteristic  of  a  singular  point  is  that  the  slope  of  the  skier  flight  path  is 
perpendicular  to  the  slope  of  the  landing  surface.  This  relationship  is  seen  in  Equation 
3.26. 

1 

tan  6  = -  (3.26) 

tan  (p 
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This  expression  can  be  written  another  way: 


(tan6)(tan(p)  +  1  =  0 


(3.27) 


Using  the  known  expressions  for  tan  9  (Equation  2.20)  and  tan  (p  (Equation 
2.25),  Equation  3.27  can  be  rewritten  into  Equation  3.28. 

/ 


tan 


tan  1  (tan  Oq  - 

\  (VoCosOoy/ 


V 


/ 


+  sin 


-1 


d^sp 


\\ 


2gh 


(»„C05  e„)2  +  iv,sin  e„  - 


(3.28) 


( tan  9o 


n 


(VoCOS  9o) 


0  +  1  =  0 


If  the  conjectured  singular  point  is  correct,  Equation  3.28  will  be  satisfied  for  any 
desired  values  of  h,Vo,  and  6q.  If  Equation  3.28  can  be  solved  with  conjecture 
coordinates  of  the  singular  point,  a  check  can  be  made  using  Equation  3.29,  which  is  the 
same  as  Equation  3.18,  but  for  a  specific  singular  point  conjecture  (Xsp,ysp)-  This 
equation  should  also  be  true  for  the  same  values  of  h  and  9q.  Table  2  summarizes  some 
results  of  this  test. 

Xsp^  —  4  cos^  00  +  4  cos^  “  4  cos^  h  tan  6q  x^p 

+  4h  cos^  00  Jsp  =  0  (3.29) 


45 


00  (degrees) 

10 

20 

30 

40 

1 

EEH  (meters) 

1 

1 

1 

1 

2 

Xsp  (conjectured) 

0.235 

0.485 

0.768 

1.114 

3 

jsp  (conjectured) 

-0.986 

-0.940 

-0.847 

-0.670 

4 

Solution  to 

Equation  3.28 

5.5511e-016 

2.2204e-016 

4.6800e-008 

1.1102e-016 

5 

Solution  to 

Equation  3.29 

9.7145e-017 

8.3267e-017 

5.5511e-016 

-2.2204e-016 

Table  2.  Test  results  to  verify  singular  point  position  for  constant  h  and  increasing  Oq 


In  Table  2,  line  4  lists  solutions  to  Equation  3.28,  which  verifies  that  the  slope  of 
the  skier  flight  path  is  perpendicular  to  the  slope  of  the  landing  surface.  Line  5  takes  the 
same  coordinates  of  the  conjectured  singular  point  and  substitutes  them  into  Equation 
3.29.  The  results,  although  not  exactly  zero,  are  so  extremely  small  that  they  can  be 
regarded  as  solutions  to  these  equations.  As  a  result,  one  can  reason  that  the  singular 
point  conjectures  are  correct. 

Eigure  31  shows  the  locus  of  the  singular  point  while  increasing  both  6q  and  h. 
As  00  increases  for  an  individual  h,  the  singular  point  moves  up  and  to  the  right  in  a 
curved  manner.  Eor  an  individual  0o  and  increasing  h,  the  singular  point  moves  down 
and  to  the  right  linearly. 
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Figure  31.  Singular  point  locations  for  /z=0.5m,  1.0m,  1.5m,  2.0m  and  0o=lO°»  20°, 

30°,  40° 

Figure  32  shows  that  keeping  6q  constant  and  increasing  h  produces  increasingly 
larger  ellipses.  These  ellipses  all  share  the  origin  as  a  common  solution.  The  singular 
point  gets  pushes  farther  to  the  right  and  down  as  h  increases.  This  makes  physical  sense, 
as  an  increased  accepted  fall  height  will  produce  a  greater  distance  from  the  skier  takeoff 
point  (origin)  to  the  singular  point. 
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Figure  32.  Oq=25  degrees  for  h=0.5m  (smallest  ellipse),  1.0m  (middle),  1.5m  (largest 

ellipse) 


C.  NON  LIPSCHITZ  CONTINUITY  OF  SAFE  SURFACE  ODE 

There  are  an  infinite  number  of  safe  landing  slope  solutions,  all  of  which  begin  at 
the  singular  point.  Notice  how,  along  with  the  uniqueness  of  the  convergence  point, 
(which  is  the  singular  point),  the  slope  of  the  bundle  of  the  safe  landing  surfaces 
asymptotically  approach  a  distinct  angle  closer  to  the  singular  point.  This  angle  is  a, 
since  at  a  specific  singular  point  (given  a  value  for  h  and  Oq),  there  exists  an  ellipse,  of 
which  the  major  axis  is  rotated  at  angle  a. 
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Figure  33.  All  safe  landing  slopes  converge  to  the  singular  point 

There  are  multiple  solutions  to  the  ODE  (Equation  2.27)  because  the  right  hand 
side  is  not  Eipschitz  continuous  at  the  singular  point. 


A  function /(t,j)  satisfies  a  Eipschitz  condition  in  the  variable  y  on  a  set  D  c 
if  a  constant  E>0  exists  with  |/(t,3/i)  — /(t, 3/2)!  <  L|yi  —  y2|,  whenever  (t,yi), 
(t,y2)  G  D.  The  constant  E  is  called  a  Eipschitz  constant  for  f  [19]. 

Close  to  the  singular  point,  there  is  a  lack  of  Eipschitz  continuity.  With  no 

Eipschitz  continuity,  there  may  a  lack  of  uniqueness  in  the  solution.  Eurthermore, 

solving  the  ODE  at  the  singular  point  can  also  lead  to  large  errors  in  finding  solutions. 

Typical  error  bounds  for  a  numerical  solution  to  an  ODE  have  a  factor  where  L 
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is  the  Lipschitz  constant  and  —  3/2  is  the  change  in  the  independent  variable.  At  the 
singular  point,  L  is  infinite,  so  the  error  can  be  huge. 

The  inverse  sine  on  the  right  hand  side  causes  the  ODE  to  not  be  Lipschitz 
continuous  at  the  singular  point.  Figure  34  shows  the  function  y  =  sin“^  x  for  (— 1  < 
X  <  1).  The  function  has  a  derivative  at  every  point  except  at  x  =  1  andx  =  —  1,  as 

■^sin~^  X  =  =  00  at  these  points.  Since  the  slope  is  infinite  at  x  =  1  and  x  =  —  1, 

the  ODE  becomes  unsolvable  when  the  inverse  sine  portion  is  1 . 


Figure  34.  y  =  sin  ^  x  for  (— 1  <  x  <  1) 

Moving  to  the  left  on  the  curve  y  =  sin“^  x  at  x  =  1  by  an  infinitesimal  amount 
makes  the  ODE  solvable. 
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IV.  IN-RUN  DESIGN 


A.  MOTIVATION  FOR  A  TRANSITION  DESIGN 

The  previous  sections  determined  a  procedure  to  develop  a  safe  landing  surface 
for  a  skier  leaving  a  jump  with  a  given  takeoff  angle.  This  section  will  outline  a  method 
to  ensure  a  safe  in-run  and  transition  portion  of  the  jump.  The  snow  surface  prior  to  the 
skier  leaving  the  jump  takeoff  can  be  modeled  with  three  distinct  parts: 

1 .  The  in-run,  a  linear  surface  that  a  skier  rides  down  from  a  given  starting 

point. 

2.  The  curved  transition  region  that  moves  the  skier  from  the  linear  surface  to 
the  linear  takeoff  ramp. 

3.  The  straight  takeoff  ramp. 

The  goal  for  this  section  is  to  determine  a  good  way  to  limit  the  rate  of  change  of 
acceleration  (jerk)  upon  the  skier,  which  can  be  accomplished  through  design  of  the 
transition  region,  based  on  a  unique  curve  known  as  a  clothoid. 

1.  Problems  With  a  Circular  Transition 

As  a  skier  moves  in  the  transition  region,  his  velocity  direction  changes  as  he 
moves  from  downhill  at  angle  —A  to  uphill  at  angle  p.  Figure  35  shows  that  a  skier  turns 
a  total  angle  of  A  +  p  prior  to  jumping. 
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In  a  circular  transition,  he  will  experience  centripetal  acceleration  with 

magnitude  — ,  r  being  the  radius  of  the  circular  path.  The  linear  in-run  surface  can  be 

considered  a  circle  with  infinite  radius,  thereby  producing  no  additional  centripetal 
acceleration.  The  normal  acceleration  is  not  continuous  at  the  point  that  a  skier  moves 
from  this  linear  slope  to  a  circular  transition.  As  a  result,  he  will  feel  infinite  jerk,  the 
time  rate  of  change  of  acceleration.  Making  the  radius  of  the  transition  circle  larger 
allows  for  the  magnitude  of  the  normal  acceleration  to  decrease,  but  even  large  radii  (50 
meters  or  more)  still  have  the  skier  feeling  infinite  jerk  at  the  two  transition  points.  This 
jerk  can  cause  balance  problems  for  the  skier  in  the  transition  and  will  not  be 
advantageous  to  the  skier’s  body  and  mind  while  preparing  to  jump  and  land  safely. 

2.  Problems  With  a  Circular  Take-Off  Ramp 

The  skier  exits  the  transition  to  a  takeoff  ramp,  which  should  be  linear  with 
angle  p.  If  a  ramp  has  any  curvature  to  it,  it  may  cause  the  skier  to  be  propelled  into  an 
undesired  back  flip  [11].  Figure  36  shows  a  terrain  park  jump  with  a  curved  takeoff 
ramp.  Consider  the  mishaps  a  skier  might  have  upon  leaving  the  jump  and  flying  through 
the  air.  The  worst  case  is  that  he  may  accidentally  rotate  backwards,  and  find  himself  in 
a  back  flip.  This  can  lead  to  him  landing  on  his  head,  which  is  undoubtedly  the  worst 
way  to  land  after  leaving  a  jump.  While  some  skiers  desire  to  complete  a  back  flip  and 
land  on  their  feet,  most  prefer  to  glide  off  of  the  jump  and  perform  some  other  type  of 
trick.  Only  a  linear  take  off  ramp  should  be  considered  for  the  jump. 


52 


Certain  types  of  ski  jumps  purposely  have  a  circular  arc  in  the  ski  jump  take  off. 
As  Figure  37  shows,  aerial  ski  jumpers  desire  a  back  flip  rotation,  as  this  is  often  a 
standard  beginning  trick  followed  by  numerous  additional  rotations  and  flips.  The 
takeoff  ramp  in  these  jumps  is  specifically  designed  to  have  curvature.  This  requires  the 
skier  to  sustain  centripetal  acceleration  upon  takeoff,  and  engenders  backward  angular 
velocity.  However,  for  typical  terrain  park  jumps  found  at  ski  areas,  most  skiers  do  not 
want  backwards  rotations  after  jumping. 


2 


Side  view  photograph  of  a  terrain  park  takeoff  ramp.  Jumps  like  this  one  exist  in  many  terrain  parks. 
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/ 

Figure  37.  A  course  layout  for  an  aerial  ski  jump,  from  [20] 

B.  DEVELOPMENT  OF  A  MODEL  FOR  A  CIRCULAR  TRANSITION 
1.  Free  Body  Diagram 

To  find  the  velocity  and  normal  acceleration  at  any  point  along  the  circular 
transition,  we  can  find  the  forces  acting  on  a  particle  (skier)  at  any  point  during  the  three 
distinct  portions  of  a  jump  prior  to  takeoff:  the  in-run,  the  curved  transition,  and  the 
linear  takeoff  ramp.  During  any  part  of  the  in-run,  the  velocity  can  be  found  through 
applications  of  Newton’s  second  law  (F  =  ma).  Since  acceleration  is  the  time  derivative 
of  velocity,  it  is  possible  to  solve  for  v  numerically  with  MATLAB. 
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Figure  38  shows  the  forees  aeting  on  a  skier  along  a  linear  in-run.  Summing 
forees  in  the  tangential  direetion  results  in  an  expression  involving  the  eoeffieient  of 
frietion  fi,  normal  foree  N,  mass  m,  gravity  g,  angle  A,  acceleration  i>,  and  variables  in  the 
drag  force:  density  of  air  p,  cross  sectional  area  A^s,  and  drag  coefficient  Q. 


mv 


—pN  +  mg  sin  A  —  -pA^sV^Cd  =  mv 


(4.1) 

(4.2) 


A  similar  expression  results  in  summing  the  forces  in  the  normal  direction,  which 
is  algebraic  since  the  normal  velocity  and  acceleration  are  both  zero. 


(4.3) 


N  —  mg  cos  A  =  0 


(4.4) 


Solving  for  N  in  Equation  4.4  and  substituting  into  Equation  4.2  results  in  a 
differential  equation  that  governs  the  velocity  of  a  skier  in  the  linear  in-run. 

1 

v  =  —pgcosA  +  gsinA  —  - — pA^sV^Cii  (4.5) 

2m 
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Integrating  Equation  4.5  with  respect  to  time  gives  the  velocity  of  the  skier  along 
the  linear  path  (Equation  4.6),  and  integrating  again  (Equation  A.lf  gives  the  position  v 
of  the  skier  at  time  t  [21]. 


V  = 


2m(g  sin  X  —  gg  cos  X) 


P-^cs^d 


tank 

f 

t 

1  ' 

pAcsCdig  sin  A  —  gg  cos  A) 


2m 


(4.6) 


2m 

X  = - In  cosh  I  t 

pAcs^d 


pA^gC^(gsinA  —  gg  cos  A) 


2m 


(4.7) 


Solving  for  t  in  Equation  4.7  and  substituting  that  in  Equation  4.6  results  in  an  expression 
for  velocity  at  the  end  of  linear  path  x  (Equation  4.8).  If  x  =  L^,  this  equation  gives  the 
velocity  at  the  end  of  a  linear  in-run. 


V 


2m(g  sin  A  —  gg  cos  A) 


P^cs^d 


tanh^cosh  ^ 


^  2m 


(4.8) 


^  Equation  4.7  was  adapted  from  an  equation  in  the  article  “Identification  of  Basketball  Parameters  for 
a  Simulation  Model”  by  H.  Okubo  and  M.  Hubbard.  In  the  article,  vertical  displacement  of  a  basketball 
from  the  floor  is  expressed  as  a  similar  expression  to  Equation  4.7.  After  applying  the  equation  in  the 
article  to  the  linear  in-run  model,  it  was  possible  to  find  the  time  derivative  solved  for  velocity  (Equation 
4.6). 
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A  similar  free  body  diagram  approach  is  used  for  the  circular  transition  (Figure  39). 


Figure  39.  Free  body  diagram  of  a  particle  in  the  circular  transition  region 


The  sum  of  forces  in  the  tangential  direction  is  equal  to  the  product  of  the  mass 
and  acceleration. 

F  =  mv 

T 

1 

—mg  sin  a  —  gN  —  -  pA^sV  =  mv 


(4.9) 

(4.10) 


Equation  4.10  uses  a  new  variable  a,  which  is  the  varying  angle  of  the  surface 
tangent  relative  to  horizontal.  The  value  of  a  ranges  from  —A  to  /S,  as  these  are  the  angles 
of  the  in-run  slope  and  take  off  ramp,  respectively. 

Similarly,  the  sum  of  forces  in  the  normal  direction  is  equal  to  the  product  of  the 
mass  and  acceleration.  In  the  circular  transition,  the  centripetal  acceleration  is  — , 
whereas  in  a  linear  portion,  the  centripetal  acceleration  is  zero. 
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(4.11) 


-mg  cos  a  -\-N  =  m- 


V 


(4.12) 


Solving  for  N  in  Equation  4.12  and  substituting  it  in  Equation  4.10  results  in 
Equation  4.13,  the  acceleration  of  a  skier  in  a  circular  transition. 

1 

V  =  —g  sin  a  —  g.( - h  g  cos  a)  -- — pAf-sV^C^  (4.13) 

r  2m 

Because  the  skier  path  is  a  circle,  velocity  is  the  product  of  the  radius  and  angular 
velocity  co. 


V  =  ro)  =  ra. 


(4.14) 


Solving  for  a  gives  the  rate  of  change  of  the  interior  and  slope  angle  a. 

V 

a=-  (4.15) 

r 

The  final  part  of  the  transition  is  the  take  off  portion.  This  free  body  diagram 
(Eigure  40)  is  similar  to  the  first  linear  section  the  in-run. 


Eigure  40.  Eree  body  diagram  of  a  particle  along  the  linear  take  off  ramp 
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Summing  forces  in  the  tangential  direction  results  in  an  expression  involving  the 
takeoff  angle  p. 

1 

—mg  sin  p  —  HN  — -pAcsV^C^  =  mVf-  (4.16) 

Summing  forces  in  the  normal  direction  gives  N. 

N  —  mgcosp=0  (4.17) 

Solving  for  N  and  substituting  into  Equation  4.16  results  in  Equation  4.18,  an 
ODE  which  governs  the  velocity  of  a  particle  along  the  linear  takeoff  ramp. 

1 

V  = —gg  cos  p  —  g  sin p  —  ^^pAcsV^Cii  (4.18) 

Just  as  Equations  4.6,  4.7,  and  4.8  show  analytic  solutions  for  velocity  along  the 
in-run  slope,  there  are  similar  analytic  solutions  for  velocity  and  position  along  the  linear 
takeoff  ramp.  However,  since  the  paths  are  not  simply  mirror  images  of  one  another,  due 
to  friction  and  drag,  it  is  easier  to  integrate  the  ODE  for  the  linear  takeoff  ramp  in  Table 
3.  Using  MATEAB  and  one  of  its  built-in  integration  functions  will  allow  us  to  find  the 
velocity  at  any  desired  point  along  the  in-run.  This  velocity  will  then  become  Vq,  the 
takeoff  velocity.  Table  3  presents  a  summary  of  equations  that  govern  velocity  for  the 
linear  in-run,  circular  transition,  and  linear  takeoff  ramp. 
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Jump  portion 

Governing  Equations 

Linear  in-run 

1 

V  =  —[ig  cos  A  +  g  sin  A  —  ^ —  pA^sV^  C^i 

which  has  solution  for  in-run  length  Lj^: 

^  ^  smA  -  ^  cos  2)^ 

Circular  Transition 

1 

V  =  —g  sin  a  —  g( - 1-  g  cos  a)  -- — pA^sV 

r  2m 

V 

a  =  — 
r 

Linear  takeoff  ramp 

1 

V  =  -pg  cosp  -  gsinp  pA^sV^ Q 

Table  3.  Differential  equations  used  to  find  the  velocity  at  any  point  in  an  in-run 


2.  Analysis  of  Circular  In-Run  Design 

Figure  41  shows  an  example  plot  of  velocity  vs.  time  in  a  circular  acceleration 
for  Q  =  0,  and  for  values  of  A=35  degrees  (incoming  slope),  P=10  degrees  (takeoff 
ramp),  /x  =  0.05  (coefficient  of  friction  between  skis  and  the  snow),  radius=7.1m, 
Lji  =  10m,  and  a  velocity  upon  entering  the  transition  of  10.22  m/s  (22.86  mph). 
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Figure  41.  Velocity  vs.  time  in  a  circular  transition,  with  A  =  35°,  =  10°,  |U=0.05, 

and  Q  =  0 


The  skier  increases  in  velocity  to  approximately  10.95  m/s  before  slowing  down. 
In  total,  he  will  spend  about  0.52  seconds  in  the  transition.  Figure  42  shows  normal 
acceleration  in  g’s  vs.  time  in  a  circular  transition. 
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Figure  42.  Ratio  of  normal  acceleration  vs.  time  in  a  circular  transition,  with  1  =  35°, 

P  =  10°,  jU=0.05,  and  Q  =  0 

In  a  circular  transition,  as  a  skier  moves  from  the  linear  slope  to  the  curve  with 
constant  radius,  he  will  experience  an  instantaneous  centripetal  acceleration.  This  jerk 
occurs  at  the  endpoints  of  the  curve  in  Figure  42.  At  the  left  endpoint,  where  the  skier 
enters  the  circular  transition,  he  experiences  an  abrupt  change  in  apparent  gravitational 
acceleration,  from  Ig  to  about  2.32g.  This  is  similar  to  having  someone  of  equal  body 
weight  suddenly  jump  on  the  skier’s  back  and  remain  there  to  the  end  of  the  transition. 
At  the  right  endpoint,  the  sudden  shift  from  2.6g  back  to  Ig  is  even  more  profound.  Keep 
in  mind  the  skier  must  now  negotiate  the  takeoff  ramp  that  is  just  a  few  feet  ahead  of  him. 
This  is  hardly  an  ideal  situation.  The  next  section  presents  a  curved  shape  that  will  limit 
the  jerk  that  a  skier  feels  at  the  transition  points,  and  will  ultimately  lead  to  a  safer 
transition  prior  to  jumping. 
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C.  DEVELOPMENT  OF  A  MODEL  FOR  A  CLOTHOIDAL  TRANSITION 

1.  The  Clothoid 

A  clothoid  is  a  unique  design  that  will  ease  the  transition  from  a  linear  slope  to  a 
curved  section.  It  has  a  variable  radius  along  its  length,  and  it  provides  a  smooth  link 
between  a  straight  line  and  a  circular  curve  [22].  It  has  the  unique  property  that  the 
distance  from  the  start  of  a  transition  s  and  the  radius  r  at  that  point  are  inversely 
proportional  to  each  other  [23].  That  is,  =  rs,  where  is  a  constant,  to  be  discussed 
later.  This  shape  at  s  =  0  initially  has  an  infinite  radius.  Clothoids  are  frequently  used  in 
rollercoaster,  highway  and  railway  design,  as  they  allow  a  car  or  train  to  travel  around  a 
loop  or  curve  without  experiencing  infinite  jerk  [23]. 

Consider  a  straight  highway  that  requires  a  necessary  curve  in  a  section.  If  the 
highway  transitions  suddenly  from  a  straight  section  to  a  circular  section,  the  vehicle  will 
experience  an  infinite  jerk  the  instant  that  it  enters  and  exits  the  circular  region  [22] .  A 
better  way  to  allow  the  road  to  bend  is  through  a  clothoid,  also  known  as  a  Cornu  spiral 
or  Euler’s  Spiral  [24].  Similarly,  if  a  rollercoaster  car  transitions  from  a  straight  to  a 
circular  section,  there  is  a  definable  point  where  the  track  changes  from  being  linear  to 
circular  with  constant  curvature  [25].  This  creates  an  undesirable  effect  on  the  riders,  as 
the  centripetal  acceleration  changes  from  zero  to  a  nonzero  magnitude  instantly.  A 
clothoid  allows  the  acceleration  to  increase  slowly  while  minimizing  stress  on  the  riders, 
yet  still  allows  the  roller  coaster  to  complete  a  loop.  The  centripetal  acceleration  varies  at 
a  constant  rate:  from  zero  force  at  the  linear  ends,  to  the  maximum  force  at  the  most 
curved  part  [25]. 

With  two  clothoids  joined  together  at  a  certain  point,  so  that  they  share  a  common 
radius,  it  is  possible  to  design  a  proper  transition  for  a  ski  jump.  Instead  of  having  the  in 
run  surface  and  takeoff  ramp  joined  by  a  circular  segment  (which  leads  to  infinite  jerk),  a 
clothoidal  transition  will  have  two  clothoids  joined  together.  One  clothoid  will  start  with 
zero  curvature  (infinite  radius)  and  turn  to  having  a  minimum  radius.  The  second 
clothoid  begins  where  the  first  one  ends,  with  minimum  radius,  and  turns  to  having  zero 


63 


curvature  (infinite  radius)  at  the  straight  slope  of  the  takeoff  ramp.  The  two  clothoids 
will  allow  the  skier  to  cope  better  with  changing  accelerations  in  the  transition  region. 

A  clothoid  can  be  defined  by  the  parametric  equations  [24] . 

m  =  f 

Jo 

=  [ 

Jo 

2 

1.5 

1 

0.5 

0 

-0.5 

-1 

-1  -0.5  0  0.5  1  1.5  2 

Figure  43.  MATLAB  plot  of  a  clothoid  (Euler’s  Spiral) 


I  [  [  c  r 


_ c _ c _ c _ t 


sin—  dx 


cos  —  dx 


(4.19) 

(4.20) 


^  The  MATLAB  code  used  in  this  plot  of  a  clothoid  was  taken  from  an  open  source  online,  found 
at  www.mathworks.in/matlabcentral/newsreader/view  thread/170988 
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The  clothoid  possesses  the  characteristic  that  the  curvature  at  any  point  is 
proportional  to  its  arc  length,  measured  from  the  origin.  In  Figure  43,  notice  that  the 
spiral  rapidly  tightens  as  the  arc  length  increases.  The  next  sections  will  discuss  the 
procedure  to  design  a  clothoid  by  using  a  tolerable  gravitational  acceleration  and  a 
minimum  transition  radius. 


2.  Minimum  Transition  Radius  Limits  Normal  Acceleration 

At  sea  level,  objects  at  rest  on  the  earth  experience  one  unit  of  the  acceleration  of 
gravity  g.  However,  humans  daily  experience  acceleration  that  are  both  greater  and  less 
than  the  force  of  gravity.  Let  y  be  the  centripetal  acceleration  in  g’s. 

The  level  of  y  can  be  arbitrarily  prescribed.  It  is  the  comfortable  level  tolerated 
by  the  human  body.  Although  humans  can  endure  forces  of  gravity  beyond  6g,  such  as  in 
a  roller  coaster  [26],  something  closer  to  1.5g  is  probably  a  much  more  reasonable 
estimate.  Although  some  skiers  can  handle  higher  sustained  gravitational  acceleration, 
this  discussion  will  focus  on  one  that  a  child  or  beginner  will  be  comfortable  with. 

The  value  of  y  can  be  thought  of  as  the  allowable  centripetal  acceleration  in  units 
of  g.  It  is  possible  to  solve  for  the  minimum  radius  of  a  curve  that  corresponds  to  y 
at  a  given  velocity.  Equations  4.21  to  4.23  present  the  steps  to  solve  for  the  radius  of 
such  a  curve. 


V 


centripetal  acceleration  = - =  gy 


V 


y  = 


^mind 

y2 

yg 


(4.21) 

(4.22) 

(4.23) 


For  any  given  desired  y  and  skier  velocity  at  the  beginning  of  the  circular 
transition,  it  is  possible  to  know  the  minimum  radius  required  to  ensure  that  the  skier 
does  not  encounter  a  higher  centripetal  acceleration  than  desired. 
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velocity  (m/s)  of  skier  at  beginning  of  transition 


Figure  44.  Circular  transition  radius  vs.  velocity  for  three  values  of  centripetal 

acceleration  y 

Figure  44  shows  the  minimum  radius  required  to  keep  a  skier  feeling  an 
additional  acceleration  of  a  desired  y  level.  Establishing  a  specific  y  allows  calculation 
of  a  minimum  transition  radius,  which  will  then  be  used  in  designing  an  exact  clothoid. 

The  arc  length  of  this  circular  transition  can  be  found  by  multiplying  the 
minimum  radius  by  the  total  angle  (A  +  /?)  that  the  skier  rotates  in  the  transition. 

5  =  r(A  +  p)  (4.24) 

Furthermore,  the  minimum  radius  required  corresponds  to  a  maximum 
approximation  for  the  transition  length  since  realistically,  the  sum  of  A  and  yS  will 
unlikely  be  greater  than  60  degrees.  (60  degrees  is  just  more  than  one  radian  =  57.29 
degrees).  As  Equation  4.24  indicates,  they  are  related  linearly. 

The  circular  length  is  found  by  Equation  4.25. 
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(4.25) 


s,=—a  +  P) 

Y9 

The  circular  length  gives  a  good  approximation  for  a  lower  bound  of  how  long  the 
clothoid  must  be.  Trials  will  show  that  a  clothoid  transition  is  approximately  twice  the 
length  of  a  circular  transition. 

3.  Development  of  Equations  to  Model  a  Clothoidal  Transition 

a.  Relationship  That  Governs  a  Clothoid 

The  relationship  that  governs  a  clothoid  can  be  expressed  as  the  product  of 
the  arc  distance  along  a  clothoid  and  the  radius  of  curvature  at  the  point  [23]. 

=  sr  (4.26) 

Where 

A  =  Flatness  of  the  curve;  a  constant 

s  =  Distance  from  the  start  of  the  transition,  where  curvature  is  zero 

r  =  Radius  of  curvature  at  distance  s  along  the  clothoidal  transition 

Figure  45.  A  family  of  clothoids  is  formed  by  fixing  the  length  s  while  varying  A  and 

r,  from  [27] 

Figure  45  shows  a  family  of  clothoids  with  varying  curvature.  In  a 
clothoidal  ski  jump  transition  design,  two  clothoids  will  join  in  the  middle  of  the 
transition  at  their  shared  point  of  maximum  curvature  and  minimum  radius.  They  are  the 
same  shape  and  the  same  length,  but  can  be  thought  of  as  offset  mirror  images  of  one 
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another.  One  clothoid  will  be  used  to  design  the  transition  from  the  in  run  surface  —A  to 
the  point  of  full  maximum  curvature  at  it’s  end.  The  second  clothoid  will  start  at  the  end 
of  the  first,  where  the  radius  is  at  a  minimum,  and  take  the  skier  to  a  desired  linear  jump 
takeoff  slope  p  with  zero  curvature.  As  a  skier  enters  the  clothoid  region  from  a  constant 
slope  —A,  there  is  a  roughly  linearly  increasing  feeling  of  centripetal  acceleration  until  he 
reaches  the  end  of  the  first  clothoid. 


b.  Clothoid  Parameters:  Turning  Angle,  Spiral  Flatness,  Clothoid 
Length,  and  Angle  From  Horizontal 

During  a  skier’s  time  in  a  clothoidal  transition,  he  will  ultimately  turn 
from  -  A  to/?,  just  as  he  would  in  a  circular  transition.  Let  ^  be  the  turning  angle  at  any 
point  along  the  clothoid.  The  differential  change  in  arc  length  is  the  product  of  the  radius 
and  the  differential  change  of  its  turning  angle. 

ds  =  r  d(  (4.27) 

For  a  clothoid,  the  constant  radius  circular  arc  is  replaced  with  a  curve  of 

varying  radius. 

dC  =  —  (4.28) 

r 

Since  =  rs,  and 

ds  s  ds 

^  ^  (4-29) 

s 

Integrating  both  sides  yields  the  turning  angle  in  terms  of  the  distance 
along  the  clothoid  and  the  spiral  flatness,  or  as  the  distance  along  the  clothoid  and  the 
associated  radius  at  that  distance. 


(  = 


2A2 


s 

2r 


(4.30) 
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With  a  given  y,  A,  and  it  is  possible  to  determine  A,  the  spiral  flatness. 
Since  (  must  turn  through  a  total  angle  A  +  ^,  the  turning  angle  of  one  of  the  (equal  and 

symmetric)  clothoids  is 


s 


(A  +  /?) 


s 


(4.31) 

(4.32) 


It  is  now  possible  to  relate  the  spiral  flatness  A  to  the  total  turning  angle 
and  length  of  each  clothoid  segment  L^. 


yl2 


i  2 

(aT^ 


(4.33) 


Since  it  is  also  possible  to  express  A^  in  terms  of  the 

minimum  radius. 


(4.34) 


We  now  use  Equation  4.33  to  solve  for  Lc,  the  total  length  of  each 
clothoid  segment. 


L/=A^(A  +  /3)  (4.35) 

L,  =  A,/(ATJj  (4.36) 


The  turning  angle  of  a  clothoid  is  the  angle  between  the  tangent  to  the 
curve  at  the  beginning  of  the  transition  and  the  tangent  to  the  curve  at  point  along  the 
clothoid.  As  shown  above,  this  turning  angle  varies  quadratically  with  the  length  of  the 
clothoid.  For  a  general  case  transition,  where  the  incoming  slope  and  jump  take  off  ramp 
do  not  have  the  same  angle,  at  the  beginning  of  the  first  clothoid  the  turning  angle  will  be 

zero.  At  the  end  of  this  clothoid,  where  5  is  at  a  maximum,  the  turning  angle  will  be 

The  turning  angle  is  the  product  of  a  constant  C  and  the  square  of  L^,  the  maximum 
length  of  the  clothoid. 
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(4.37) 


A  +  p 
2 


=  CL,^ 


After  solving  for  C  in  Equation  4.37,  the  turning  angle  can  now  be  found 
at  any  value  of  s  along  the  clothoid. 


^  =  Cs^ 


(A  +  ^)  2 

- 


(4.38) 


In  the  first  clothoid,  a  skier  enters  the  transition  at  a  point  of  infinite  radius 
(zero  curvature)  and  travels  to  a  point  of  known  minimum  radius,  at  the  end  of  the  first 
clothoid.  As  he  gets  closer  to  the  end  of  this  clothoid,  his  rate  of  turning  increases.  At 
the  end  of  the  first  clothoid,  he  enters  the  second  clothoid.  This  clothoid  begins  at  the 
point  of  minimum  radius,  and  ends  at  a  point  of  infinite  radius,  which  is  the  beginning  of 
the  takeoff  ramp  with  constant  slope  /?.  As  a  result,  the  skier  feels  a  gradual  centripetal 
acceleration,  as  opposed  to  in  a  circular  transition,  where  there  is  an  instant  of  infinite 
jerk. 

Another  meanginful  characterization  of  a  clothoid  transition  is  based  on  0, 
the  angle  of  the  surface  from  horizontal.  In  the  transition,  at  the  point  of  intersection 
between  the  linear  in-run  surface  and  the  clothoid,  =  0  and  0  =  —A.  At  the  intersection 

of  the  two  clothoids  d  =  and  0  =  At  the  end  of  the  second  clothoid,  where  the 
^22 

entire  transition  ends,  the  value  of  (  =  A  +  p,  and  the  value  of  0  =  yS.  A  summary  of 
these  angles  is  seen  in  Figure  46  and  Table  4. 

The  average  rates  of  change  of  (  and  0  are  the  same. 


dS  =  (4.39) 

Integrating  both  sides  yields  an  expression  relating  the  turning  angle  and 
the  angle  from  horizontal. 


0  =  <■  +  C  (4.40) 

The  value  of  C  (constant  of  integration)  is  found  from  the  initial  condition 
of  0  =  —A  when  ^  =  0.  Therefore,  everywhere  along  both  clothoid  segments,  0  =  ^  — 


70 


A.  In  general,  the  angle  from  horizontal  at  any  point  along  a  clothoid  equals  the  total 
turning  angle  at  any  point  plus  the  originating  linear  surface  angle. 


Figure  46.  Relationship  between  turning  angles  and  angles  from  horizontal 


Clothoid  1 

Clothoid  2 

Initial 

Final 

Initial 

Final 

Turning  Angle  ^ 

0 

x  +  p 

2 

x  +  p 

2 

x  +  p 

Angle  from  horizontal  0 

-A 

X-f 

2 

X-(3 

2 

Table  4.  Summary  of  turning  angles  and  angles  from  horizontal  in  a  clothoid  transition 


c.  Free  Body  Diagram  of  a  Skier  in  Clothoidal  Transition 

In  a  clothoidal  transition  region,  the  free  body  diagrams  for  the  two  linear 
portions  (Figures  38  and  40)  are  the  same  as  in  a  circular  transition.  The  distinction  is  in 
the  transition,  which  is  now  clothoidal.  Figure  47  shows  the  free  body  diagram  for  a 
clothoidal  transition. 
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Minimum  radius 


mg 


Figure  47.  Free  Body  Diagram  for  a  clothoidal  transition 


Summing  the  forces  in  the  tangential  direction  produces  an  expression 
involving  the  angle  from  horizontal. 


mv 


—mg  sin  0  —  gN  —  —pAcsV^C^  =  mv 


(4.41) 

(4.42) 


Summing  forces  in  the  normal  direction  produces  a  similar  expression. 

\  F  =  mv  (4.43) 

4-1  N 

(4.44) 


17" 

—mq  cos  0  +N  =  m  — 
r 

Solving  for  N  in  the  Equation  4.44  and  substituting  it  in  Equation  4.42 
yields  the  following  expression: 

..2 


v‘-  1 

V  =  —g  sin  0  —  g( - h  g  cos  0)  —  - —  pA^sV  Q 

r  2m 


(4.45) 
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Along  the  path  of  a  clothoid,  the  velocity  of  a  particle  is  the  rate  of  change 
of  its  position  at  that  point. 


(4.46) 


d.  Clothoidal  Transition  Velocity  Equations 

For  the  first  clothoid,  the  formula  for  the  velocity  at  any  point  in  the 
clothoid  can  be  written  as  seen  in  Equation  4.45.  This  can  be  rewritten  using  the  turning 
angle  angle  A,  position  along  the  clothoid  5,  spiral  flatness  A,  density  of  air  p,  cross 
sectional  area  A^s,  and  drag  coefficient  Q. 

V  =  A)  +  g  cos{?,  -  X))  -^pA^r^Q  (4.47) 

And  finally,  using  only  v,  s.  A,  A,  p,  A^.^,  and  Q. 

V  =  -g  sin(^  -  X)  -  p(^v^  +  g  cos(^  -  A))  -  —pA.sV^C^  (4.48) 

For  the  second  clothoid,  if  the  turning  angle  ^  is  measured  from  the 
beginning  of  the  second  clothoid.  Equation  4.45  can  also  be  rewritten  using  turning 
angle  takeoff  ramp  angle  p,  position  along  the  clothoid  5,  spiral  flatness  A,  density  of 
air  p,  cross  sectional  area  A^s,  and  drag  coefficient  Q. 

v  =  -g  sin((  +  /?)  -  +  g  cos«'  +  P))  -  ^pA„v^Q  (4-49) 

And  finally,  using  only  v,  5,  p,  A,  p,  A^s,  and  Q. 

5^  S  5^  1 

V  =  -g  sin(^  +  P)  -  p(-^v^  +  g  cos(^  +  p))  - —pA^sV^Ca  (4-50) 

The  differential  equations  for  the  velocity  (Equations  4.47  and  4.49)  and 
the  rate  of  change  of  the  turning  angle  (Equation  4.51)  make  it  possible  to  find  the 
velocity  of  a  skier  at  any  point  in  the  clothoidal  transition  region  using  an  automated 
integration  function  in  MATEAB. 
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(4.51) 


d(  =  c  = 


s  ds 


A  summary  of  equations  that  govern  velocity  in  a  clothoidal  transition  is 


seen  in  Table  5. 


Jump  portion 

Governing  Equations 

Linear  in-run 

1 

V  =  —[ig  cos  A  +  g  sin  A  —  - —  pA^sV^  Q 

which  has  solution  for  in-run  length 

(  \2m(gsinA-  ggcosA)\  /  /  ^ 

r-M  \tanh[cosh^[e-i  2n, 

Clothoidal 

Transition 

1 

V  =  —g  sin  0  —  p( - 1-  g  cos  0)  -- — pA^sV 

r  2m 

ds 

V  =  —  =  s 
dt 

s  s 

Linear 

takeoff  ramp 

1 

V  =  -pg  cosf  -  gsinf  -  —pAcsV^C^ 

Table  5.  Summary  of  velocity  ODEs  for  a  clothoidal  transition 


e.  Physical  Coordinates  of  a  Clothoid 

The  (x,y)  coordinates  of  any  point  on  the  clothoid  can  be  expressed  as 
functions  of  s,  the  distance  at  any  point  along  a  clothoid  [28].  The  (x,y)  coordinates  when 
plotted,  result  in  a  clothoid  that  starts  at  the  origin,  beginning  with  zero  curvature. 

,5 


S-'  S' 

^  ^  ^  ~  40A^  ^  3456^2  ^ 
^11 


5^  5^ 


y  = 


+ 


6212  336^6  422402110 


■  + 


(4.52) 

(4.53) 
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Neither  of  the  two  conjoined  clothoids  satisfy  the  conditions  met  in 
Equations  4.52  and  4.53.  The  {x,y)  coordinates  must  be  translated,  rotated,  and  reflected 
to  give  the  true  coordinates  that  are  needed  in  building  the  correct  clothoid  transition. 
Consider  the  following  example  clothoid  plots,  all  which  use  7=1.5  and  L;j=10  meters. 


X  (m) 

Figure  48.  Complete  transition  composed  of  two  clothoids,  with  X  =  35°,  p  =  10°, 

7=1.5,  L^=10m,  Q=0,  |U=0.05 

In  Figure  48,  the  clothoid  transition  begins  at  1  =  35°  and  ends  at  ^  =  10°. 
This  plot  is  actually  two  clothoids,  joined  together  at  the  point  with  minimum  radius, 
which  is  determined  by  skier  velocity  and  desired  7  level.  Figures  49  and  50  show  other 
examples  of  clothoid  transitions  for  different  values  of  A.  and  p. 
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4.  Design  Algorithm  for  Finding  a  Clothoidal  Shape 

A  design  algorithm  for  finding  the  needed  clothoid  shapes  can  be  summarized 
with  this  algorithm: 

1.  Specify  the  degree  of  the  incoming  slope  A  (in  degrees)  and  the  desired 
takeoff  ramp  slope  angle  p  (in  degrees). 


2.  Specify  the  length  (in  meters)  of  the  incoming  slope 

3.  Determine  the  velocity  of  the  skier  upon  entering  the  transition, 


V 


2m(g  sin  A  —  (ig  cos  A) 


tr 


pA^cs^d 


tank  ycosh 


P-^cs^d 

2m 


4.  Establish  y,  the  lateral  acceleration  in  g’s  that  is  can  be  tolerated  by  a 
skier.  Realistically,  the  value  should  not  be  greater  than  2.0. 

y  2 

5.  Determine  the  minimum  radius  required,  with 

6.  Find  A,  the  spiral  flatness,  (constant  clothoid  design  parameter),  from 
A  =  r^in^iA  +  P). 

7.  Find  the  total  length  of  each  clothoid,  by  =  A-^{A  +  P). 

8.  Draw  the  generic  (x,y)  coordinates  of  a  clothoid  found  by  using  a 
uniformly  increasing  s  value  and  A. 

^  3456^8^’" 

^  ^  ~  336A6  ^  42240A10  ^ 

9.  Translate  and  rotate  both  clothoid  segments,  and  reflect  only  the  second 
clothoid  to  calculate  the  (x,y)  coordinates  of  the  entire  clothoidal  transition  region. 
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5.  Analysis  of  a  Clothoidal  In-Run  Design 

Figure  51  shows  that  the  pattern  of  velocity  vs.  time  in  a  clothoid  transition  is 
very  similar  to  that  of  a  circular  transition.  The  velocity  increases  to  about  1.05  m/s  more 
than  in  a  circular  transition,  and  the  time  spent  in  the  clothoid  is  nearly  double  that  of  a 
circle. 


Figure  51.  Velocity  vs.  time  for  a  clothoidal  transition  for  1=35°,  ;S=10°, 

H=0.05,  y  =  1.5,  and  Q  =  0 

Figure  52  shows  that  in  a  clothoidal  transition,  the  normal  acceleration  increases 
gradually  and  then  decreases  gradually  after  about  halfway  through  the  transition.  This  is 
because  the  unique  shape  of  the  clothoid  takes  a  skier  from  the  linear  in-run  (with  around 
Ig)  to  a  point  of  minimum  radius,  where  the  gravitational  acceleration  matches  that  felt  in 
a  circular  transition,  and  then  back  to  a  linear  surface  (the  takeoff  ramp). 
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time,  (s) 


Figure  52.  Ratio  of  normalized  mass  specific  force  to  g  vs.  time  in  a  clothoid 
transition,  for  A  =  35°,  ^  =  10°,  jU=0.05,  y  =  1.5,  and  Q  =  0 


D.  ADVANTAGES  OF  A  CLOTHOIDAL  TRANSITION  OVER  A 

CIRCULAR  TRANSITION 

The  clothoid  shaped  transition  is  an  improvement  over  a  circular  transition,  as  the 
skier  experiences  a  gradual  change  in  apparent  normal  acceleration  rather  than  one  with 
sudden  changes  in  acceleration  and  infinite  jerk.  Considering  a  specific  case  of  defined 
parameters,  an  analysis  of  velocities  and  forces  in  the  transitions  can  be  conducted. 

Notice  how,  in  a  clothoidal  transition,  the  time  spent  in  the  transition  is  nearly 
double  that  for  the  circular  transition.  This  is  because  the  clothoid  shaped  transition  is 
about  twice  the  length  of  a  circular  shaped  transition.  In  both  plots,  the  skier  enters  the 
transition  from  the  linear  existing  slope  at  the  same  speed,  and  the  skier  experiences  an 
increase  in  velocity,  then  an  instant  where  the  velocity  begins  to  decrease.  The  velocity 
will  decrease  throughout  the  remainder  of  the  transition.  Although  these  velocity  plots 
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(Figures  41  and  51)  are  useful,  a  comparison  of  the  normal  forces  (Figure  42  and  52)  in 
the  transitions  will  uncover  the  advantage  that  a  clothoid  transition  has  over  a  circular 
transition. 

In  a  circular  transition,  as  a  skier  moves  from  the  in-run  to  the  beginning  of  the 
circular  arc,  he  will  encounter  infinite  jerk  (Figure  53).  From  -2  to  0  seconds,  the  skier  is 
on  the  in-run.  From  0  to  0.5  seconds,  the  skier  is  in  the  circular  transition,  and  at  0.5 
seconds,  the  skier  begins  moving  along  the  linear  takeoff  ramp.  Notice  the  infinite  jerk 
that  the  skier  experiences  at  the  ends  of  the  circular  transition,  represented  by  the  vertical 
lines. 


Figure  53.  Ratio  of  normal  acceleration  vs.  time  in  the  in-run,  circular  transition,  and 
takeoff  ramp  for  A=35°,  ^=10°,  jU=0.05,  y  =  1.5,  and  Q  =  0 
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Figure  54.  Ratio  of  normal  acceleration  vs.  time  in  the  in-run,  clothoidal  transition, 
and  takeoff  ramp  for  A=35°,  ^=10°,  n=0.05,  y  =  1-5,  and  Q  =  0 

In  a  clothoidal  transition,  as  a  skier  moves  from  the  in-run  to  the  beginning  of  the 
clothoidal  arc,  he  will  encounter  a  measureable  amount  of  jerk.  The  lateral  acceleration 
begins  just  under  Ig,  and  increases  to  just  over  3g,  where  it  decreases  again  to  the 
previous  value.  In  Figure  54,  from  -2  to  0  seconds,  the  skier  is  on  the  in-run.  From  0  to  I 
second,  the  skier  is  in  the  transition,  and  at  I  second,  the  skier  begins  moving  along  the 
linear  takeoff  ramp. 
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Figure  55.  Comparison  of  the  ratio  of  normal  acceleration  vs.  time  for  circular  and 
clothoidal  transitions  in  the  in-run,  transition,  and  takeoff  ramp  for 
X  =  35°,  ^  =  10°,  jU=0.05,  7  =  1.5,  and  Q  =  0 

Figure  55  shows  an  overlay  of  Figure  53  and  54.  In  the  circular  transition,  a  skier 
encounters  an  instantaneous  jump  of  approximately  2.5g  at  the  beginning  and  the  end  of 
the  transition.  Conversely,  in  a  clothoidal  transition,  a  skier  does  not  have  any 
instantaneous  jump  at  all.  Rather,  the  normal  acceleration  builds  to  a  reasonable  level, 
and  then  decreases  in  roughly  the  same  manner.  Following  each  transition,  (0.5  seconds 
for  the  circular  and  1  second  for  the  clothoidal)  the  skier  begins  moving  along  the  linear 
takeoff  ramp.  Realistically,  a  skier  will  not  be  on  a  takeoff  ramp  for  as  long  as  1  or  1.5 
seconds,  but  this  figure’s  purpose  is  to  show  how  the  two  transition  shapes  allow  for 
different  methods  to  handling  jerk  along  a  vertical  curve. 

A  clothoidal  transition  will  better  enable  skiers  to  cope  with  the  centripetal 
acceleration  that  must  be  endured  while  transitioning  from  a  linear  in-run  to  a  takeoff 
ramp. 
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V.  TERRAIN  PARK  JUMP  MAINTENANCE 


A.  NEED  FOR  JUMP  MAINTENANCE  PROGRAM 

In  order  to  ensure  that  these  jumps  consistently  allow  the  skier  to  jump  safely,  a 
maintenance  program  must  be  in  place.  If,  through  repeated  landings  in  the  same  spot, 
the  safe  landing  slope  develops  a  sort  of  pothole  appearance,  the  jump  should  be  closed 
while  repairs  are  made.  The  landing  area  should  be  level  across  the  width  of  its  surface, 
as  all  landings  should  have  an  even  surface.  Currently,  little  is  known  about  how 
snowfall  affects  the  surface  of  a  ski  jump  transition,  ramp,  and  landing  surface. 
However,  this  section  gives  a  short  analysis  of  snowmelt  and  fall. 


B.  SNOWMELT 


The  amount  of  snow  that  melts  can  be  determined  by  a  simple  expression  [29]. 

M  =  c(r,-r,)  (5.1) 

with 


M  =  Snowmelt,  inches  per  day 

in 

C  =  Melt  rate  coefficient  Cttt-; — ) 

^°F  day^ 

Ta  =  Air  temperature,  °F 
Tfj  =  Base  temperature  at  which  snow  melts,  °F 

The  melt  rate  factor  typically  varies  between  0.04  and  0.08  in/°F  day,  and  the  base 
temperature  is  generally  set  to  32  °F.  Figure  56  shows  a  plot  of  snowmelt  per  day  vs. 
temperature. 

This  simple  model  is  probably  adequate  for  terrain  park  snowmelt  analysis.  It 
does  not  take  any  additional  snowmelt  induced  by  wind  or  rain,  but  it  provides  a 
conservative  estimate  of  the  amount  of  snow  that  melts  by  air  temperature  alone.  The 
plot  below  gives  snowmelt  estimates  for  various  typical  values  of  snowmelt  rate 
coefficients,  as  the  ambient  air  temperature  increases.  Naturally,  for  snow  to  melt,  the 
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temperature  must  be  above  freezing.  The  higher  the  outside  temperature  is,  the  greater 
that  the  snowmelt  will  be.  It  must  be  noted  that  at  a  reasonable  outside  temperature,  such 
as  50  degrees,  the  ski  area’s  overall  snowpack  will  not  be  sustainable,  and  conditions  will 
rapidly  deteriorate.  For  a  nice  spring  temperature  of  50  degrees  °F,  with  a  high  melt  rate 
coefficient,  snow  will  melt  at  about  1.5  inches  per  day.  As  the  daily  temperatures 
generally  get  warmer,  with  this  fast  rate,  the  ski  area  might  begin  to  close  runs  and  shut 
down  winter  operations. 


Figure  56.  Snowmelt  per  day  vs.  temperature 


In  freezing  temperatures,  the  coefficient  of  friction  of  skis  against  the  snow  will 
be  around  0.05.  For  temperatures  above  freezing,  the  apparent  coefficient  of  friction  will 
increase,  as  the  snow  becomes  slushy. 
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C.  SNOWFALL 

As  new  snow  falls  and  builds  depth  on  a  terrain  park  jump,  special  care  must  be 
taken  to  ensure  the  jumps  are  maintained.  A  terrain  park  manager  should  ensure  that  the 
additional  snow  is  packed  down  equally  throughout  the  surface  of  the  transition,  takeoff 
ramp,  and  landing  surface,  so  that  the  takeoff  angle  and  safe  slope  surface  remain 
relatively  constant.  Little  is  known  about  how  heavy  snowfall  affects  the  surface 
dimensions  and  characteristics  of  a  ski  jump. 
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VI.  FINAL  ALGORITHM 


In  determining  a  jump  design,  there  is  freedom  in  defining  many  parameters. 
First,  as  there  is  an  infinite  family  of  safe  landing  surfaces,  there  is  freedom  in 
determining  how  long  and  high  the  landing  surface  should  be.  Although  there  is  a 
restriction  on  the  degree  of  the  existing  linear  in-run  surface,  there  is  then  freedom  in 
choosing  the  takeoff  angle. 

1.  Find  a  region  on  a  ski  hill  where  one  wishes  to  build  the  jump. 

2.  Determine  the  mean  slope  of  the  hill  A  and  establish  an  in-run  length 

3.  Determine  the  desired  maximum  height  and  length  of  the  entire  jump  at  the 
highest  velocity.  These  will  be  the  (Xi,  yi)  constraints. 

4.  Determine  the  desired  EFH  (ideally  should  be  0.5  to  1.5  meters) 

5.  Choose  a  reasonable  takeoff  angle  /?. 

6.  Run  the  clothoid  algorithm  for  the  chosen  angle.  Determine  the  velocity  on 

takeoff. 

7.  Calculate  a  set  of  safe  jump  landing  surfaces  for  the  chosen  takeoff  angles  that 

fit  within  the  limitations.  These  safe  landing  surfaces  all  have  the  same 

horizontal  extent  but  vary  in  height  at  their  end. 

8.  With  this  set  of  landing  surfaces,  determine  where  to  best  place  one  on  the 
chosen  region  on  the  hill.  This  can  be  accomplished  by  shifting  the  set  of  safe  landing 
surfaces  so  that  the  area  between  the  existing  slope  and  safe  landing  surface  is 
minimized.  Choose  the  one  safe  surface  that  minimizes  this  area.  This  will  result  in  the 
least  movement  of  snow  to  build  the  landing  surface. 
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APPENDIX.  RELEVANT  MATLAB  CODES 


The  MATLAB  Codes  presented  in  this  section  are  actual  codes  used  to  produce 
the  results  in  this  thesis. 

A.  SAFE  LANDING  SURFACE  CODES 
1.  Impact  Events 

% -  Function  events  - 

function  [value,  isterminal,  direction]  =  impactevents (t, y) 
global  xs  ys  dysdxs 

%  There  are  two  events:  to  stop  integrating 

%  1.  When  the  skier  hits  the  slope,  stop  integrating 

%  2.  When  the  slope  of  the  jumper  path  =  slope  of  landing  slope 

%  When  the  jumper  is  at  the  max  air  height 

%  slope  of  jumper  path=slope  of  landing  surface  =  max  height 

format  long 

%states  z  =[  X  vx  y  vy]  These  are  the  states  required  for  ODE45 
%if  t>  0.05 

xtraj=y{l);  %  distance  in  the  x  direction 
vxtraj=y (2) ;  %  velocity  in  the  x  direction 

ytraj=y{3);  %  distance  in  the  y  direction 
vytra j =y { 4 ) ;  %  velocity  in  the  y  direction 

ty=[t  y'];  %  establish  a  vector 

yslope=spline (xs , ys , xtra j ) ; 

value  (1)  =  ytraj -yslope;  %  path  height  above  surface 

%  jumper  height  -  yslope 
%  when  this  goes  to  zero,  the  event  happens. 

%  Looks  for  a  time  when  quantity  goes  to  zero 
isterminal { 1 )  =  1; 

%  when  this  happens,  it  tells  ODE45  to  stop  integrating 
direction (1)  =  -1; 

%  when  this  happens,  program  looks  for  value  of  when  solution 
%  crosses  0  to  1 
dydxtra j =vytra j /vxtraj ; 

%dydxslope=interpl (xs , dysdxs , xtra j , 'spline ' ) 
if  xtraj  >xs { 1 ) 

dydxslope=interpl (xs , dysdxs , xtra j ) ; 

else 

dydxslope=dysdxs (1) ; 

end 

value (2 ) =dydxtraj -dydxslope;  %  slope  of  path  wrt  snow  slope 
%  =  difference  between  trajectory  derivative  and  slope  derivative 
%  helps  determine  vertical  distance  where  the  skier  will  be 
isterminal (2) =0; 
direction{2)=-l; 

results_impactevents= [xs { 1 )  xtraj  dydxtraj  dydxslope  value{2)]; 
format  short 
return 
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2.  Jump  Events 


% -  Function  events  - 

function  [value,  isterminal,  direction]  =  jumpevents (t, z) 
global  Xend 

%states  z  =[  X  vx  y  vy]  These  are  the  states  required  for  ODE45 
X  =  z  { 1 ) ;  %  distance  in  the  x  direction 

%  FINISH  EVENT 

value (I)  =  Xend-x;  %  distance  to  go  to  path  end 
isterminal { 1 )  =  1; 
direction (1)  =  -1; 
return 

%  integrates  until  x  equals  end  value  that  we  want 
%  Xend  is  in  global 

3.  Initial  Velocity 

function  [vo  v  theta] =initvel (x, y) 

%Initial  velocity  of  the  skier  before  the  skier  jumps  to  reach  point 
(x,  y) 

%Includes  calculation  for  drag  and  skier  pop 
%  once  the  skier  gets  there,  what  is  x  and  y? 

%  iterate  through  while  loop 
format  compact 

global  rho  area  mass  g  del  v  betao  sbo  cbo  tbo; 
global  Xend 

xandy  initvel=[x  y]  ; 
sbo=sin (betao) ; 
cbo=cos (betao) ; 

tbo=tan (betao) ;  %  line  11  is  drag  free  velocity 
vo=sqrt (x^2*g/ (2*cbo^2* (x*tbo-y) ) )  ;  %velocity  of  jumper  just 

before  jumping 

dely=l;  %preallocates  value  of  jumper  height  to  slope 

%height 

%  at  value  of  x 

%  integrate  jump  flight;  calculate  y 
%  if  not  the  right  y,  go  back 
%  defy  from  last  time  and  the  one  I  want. 
i=0;%allows  for  increment  change  in  line  49 
while  abs (defy) >  l.Oe-5  %runs  loop  until  moment  of  impact 
i=i+l;  %  counter 

tspan= [ 0 : 0 . 1 : 5 ] ;  %betao  is  increased  by  delta_theta 
delta  theta=atan (del  v/vo) ;  %  takeoff  angle  increment  after  pop 
thetao=betao  +  delta  theta;  %total  initial  jumper  angle  with 
horizontal 

%  jumper  velocity  is  vj 

vj=(vo^2+del  v^2)^.5  ;  %calculates  the  jumper  velocity 

%adding  skier  pop 
%component 

z0=[0  vj *cos (thetao)  0  vj *sin (thetao)  ];  %  state  variable  x, 

vx,  y,  vy 

options  =  odeset (' Events ', @ jumpevents) ;  %look  for  event 
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Xend=x;  %jumpflight  has  drag  already  in  it  %x  is  value  we  want, 
with  drag 

[t, z] =ode45 {@jumpflight, tspan, zO, options) ;  %already  includes  pop 
tz=[t  z];%solves  for  vx,  the  x  component 

%of  skier  velocity  including  drag 
xi=z {:,!);  %x  y  vxi  vyi  are  columns  of  z 
yi=z { : , 3) ; 
vxi=z { : , 2 ) ; 
vyi=z  {  : , 4)  ; 
xiyi= [xi  yi ] ; 

yact  =  z{end,3);  %  actual  value  of  y  is  third  column  of  z. 
vxa  =  z (end, 2 ) ; 
vya  =  z (end, 4 ) ; 

theta=atan (vya/vxa) ;  %  at  end  value 
v=sqrt {vxa^2+vya^2 ) ;  %value  of  y  at  v 

dely=yact-y;  %  if  vo  is  not  correct,  then  y  is  not  correct. 

%yact  is  less  than  y. 
if  i==l  %  first  time  through 

d_dely_dvo=x*x*g/ {vo^3*cos (thetao) ^2 )  ;  %  derivative 

else 

del_dely=dely-delyold; 
d_dely_dvo=-del_dely/ delvo; 
end  %  change  vo .  Divide  dely  we  need  by  d  dely  dvo 
delvo= {dely/d_dely_dvo)  ; 
vo=vo-delvo; 
delyold=dely; 

end 

4.  Safe  Surface 

function  yprime=saf esurf ace (x, y)  %  output  is  yprime;  input  is  {x,y) 

%  this  routine  implements  the  differential 
%  equation  satisfied  by  a  ski  jump  surface  of 
%  constant  equivalent  fall  height  h  for  all  speeds  vo 
%  with  a  takeoff  ramp  of  tho  degrees  assuming  no  drag 
global  h  g  betao  %  tho  cto  sto  tto 
%  betao  is  ramp  angle;  not  jump  angle 
xandy_safesurf ace= [x  y] ; 
cons=[h  g  betao];  %  cto  sto  tto  tho 
sbo=sin (betao) ; 
cbo=cos (betao) ; 
tbo=tan (betao) ; 

vo=sqrt (x^2*g/ (2*cbo^2* (x*tbo-y) ) ) ;  %initial  velocity  needed  to  impact 
at  X,  y 

%  with  no  drag  or  lift 
%  same  as  eqn  (2)  on  page  4 

%  v=total  velocity.  The  perpendicular  component  is  sqrt(2gh)  for  a 
given  h 

[vor  V  theta] =initvel (x, y) ;  %  routine  to  iterate 

%  init  V  reqd  to  get  to  x  and  y 

%  v=total  velocity.  The  perpendicular  component  is  sqrt(2gh)  for  a 
given  h 

xy=[x  y] ;  %  establishes  a  vector  of  x  and  y 

%  theta  is  the  impact  angle  between  the  velocity  vector  and  horizontal 
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V 


%  same  as  eqn  (4)  on  page  4 
if  V  <  sqrt{2*g*h) 
v=sqrt {2*g*h) ; 

end; 

zeta=-asin (sqrt {2*g*h) /v) ; 

%  zeta  is  the  allowed  angle  between  slope  and  path  at  speed 

%  zeta  is  the  angle  whose  sine  is  sqrt{2gh)/v  =  v{perp)/v 

%  therefore,  sin  { zeta) =v (perp) /v 
%  if  V  is  too  small,  then  sqrt  {2gh)  is  greater  than  v,  and  imaginary 
phi= (theta-zeta) ;  %  gamma  =  theta-beta 

%  phi  is  the  angle  between  the  slope  and  horizontal  at  an  impact  point 
yprime=tan (phi ) ;  %  surface  slope  relative  to  horizontal 

%  If  the  direction  of  the  velocity  vector  is  known,  and 
%  the  magnitude  at  impact  is  known,  and 

%  the  angle  between  v  and  the  slope  is  known,  then  we  can  find  out  how 

the 

%  slope  should  be  oriented, 
yprime { 1 ) ; 

format  long; 

rssur=[x  y  vor  v  theta  zeta  yprime]  ; 
format  short; 


5.  Safe  Jump 

%  Original  code  by  Mont  Hubbard,  1  March  2007 
%  Additional/changed  code  by  Andrew  Swedberg,  24  July  2009 
%  This  program  calculates  the  shape  of  a  ski  jump  landing  surface 
%  that  results  in  a  constant  equivalent  fall  height  at  the  impact  point 
%  of  skier  trajectory  on  the  jump  surface  as  a  function  of  takeoff 
angle . 

g, 

o 

%  The  design  uses  a  drag  free  approximation  and  the 
%  concept  of  the  zenith  in  the  design  of  the  safe  landing  slope 
%  The  slope  is  tested  using  a  simulation  including  drag  however. 

g, 

o 

clear 

clf 

hold  off 

global  h  g  betao  rho  mass  area  del  v  %  tho  cto  sto  tto 

global  ft2m  m2ft  deg2rad  mps2mph  mph2mps 

global  hnom 

global  xs  ys  dysdxs 

g, 

o 

format  compact;  hold  off; 

rad2deg=180/pi;  %  converts  radians  to  degrees 
ft2m=0.3048;  %  converts  feet  to  meters 

m2ft=l/ft2m;  %  converts  meters  to  feet 

deg2rad=pi/180; %  converts  degrees  to  radians 

mps2mph=2 . 2369362 9;  %  converts  meters/second  to  miles/hour  (changed 
from  2.237) 

mph2mps=l/mps2mph; %  miles  per  hour  to  meters  per  second  conversion 
g=9. 80665;  %  meters  per  second  (changed  from  9.81) 

96 


g, 

o 

%  to  convert  lb/ft3  to  kg/mS,  multiply  by  16.0185 
%  to  convert  kg/m3  to  lb/ft3,  multiply  by  0.0624 
%  Squaw  Valley  Base  elevation:  6,200  ft  (1,890  m) 

%  Squaw  Valley  Top  elevation:  9,050  ft  (2,760  m) 

%  Aspen  Highlands  Base  elevation:  8,040  feet  (2,450  m) 

%  Aspen  Highlands  Top  elevation:  11,678  feet  (3,559  m) 

%  Jackson  Hole  Base  elevation:  6,311  feet  (1,924  m) 

%  Jackson  Hole  Top  elevation:  10,450  feet  (3,190  m) 

g, 

o 

%  From:  http://www.denysschen.com/catalogue/density.asp 

g, 

o 

%  Uncomment  any  rho  US  desired,  depending  on  elevation  of  jump 
%  rho  US  is  measured  in  lb/ft^3  ;  see  conversion 

%  rho_US=0 . 0645;  %  For  5000  feet  at  (IOC,  50F) 

%  rho_US=0 . 0658;  %  For  5000  feet  at  (4.44C,  40F) 

%  rho_US=0 . 0668;  %  For  5000  feet  at  (OC,  32F)  freezing 

%  rho_US=0 . 0685;  %  For  5000  feet  at  (-6.67C,  20F) 

%  rho_US=0.07;  %  For  5000  feet  at  (-12.22C,  lOF) 

%  rho_US=0 . 07 15 ;  %  For  5000  feet  at  (-17.77C,  OF) 

%  rho_US=0.0731;  %  For  5000  feet  at  (-23.33C,  -lOF) 

%  rho_US=0 . 0621;  %  For  6000  feet  at  (IOC,  50F) 

%  rho_US=0.0633;  %  For  6000  feet  at  (4.44C,  40F) 

%  rho_US=0 . 0644 ;  %  For  6000  feet  at  (OC,  32F)  freezing 

%  rho_US=0 . 066;  %  For  6000  feet  at  (-6.67C,  20F) 

%  rho_US=0.0674;  %  For  6000  feet  at  (-12.22C,  lOF) 

%  rho_US=0.0688;  %  For  6000  feet  at  (-17.77C,  OF) 

%  rho_US=0.0704;  %  For  6000  feet  at  (-23.33C,  -lOF) 

g, 

o 

%  rho_US=0 . 05 98 ;  %  For  7000  feet  at  (IOC,  50F) 

%  rho_US=0.061;  %  For  7000  feet  at  (4.44C,  40F) 

%  rho_US=0 . 062 ;  %  For  7000  feet  at  (OC,  32F)  freezing 

%  rho_US=0 . 0635;  %  For  7000  feet  at  (-6.67C,  20F) 

%  rho_US=0.0649;  %  For  7000  feet  at  (-12.22C,  lOF) 

%  rho_US=0 . 0663;  %  For  7000  feet  at  (-17.77C,  OF) 

%  rho_US=0.0678;  %  For  7000  feet  at  (-23.33C,  -lOF) 

g, 

o 

%  rho_US=0 . 0575 ;  %  For  8000  feet  at  (IOC,  50F) 

%  rho_US=0.0587;  %  For  8000  feet  at  (4.44C,  40F) 

rho_US=0 . 05 97 ;  %  For  8000  feet  at  (OC,  32F)  freezing 

%  rho_US=0 . 0611 ;  %  For  8000  feet  at  (-6.67C,  20F) 

%  rho_US=0.0624;  %  For  8000  feet  at  (-12.22C,  lOF) 

%  rho_US=0.0638;  %  For  8000  feet  at  (-17.77C,  OF) 

%  rho_US=0 . 0652 ;  %  For  8000  feet  at  (-23.33C,  -lOF) 

%  rho_US=0 . 0554 ;  %  For  9000  feet  at  (IOC,  50F) 

%  rho_US=0.0565;  %  For  9000  feet  at  (4.44C,  40F) 

%  rho_US=0 . 0574 ;  %  For  9000  feet  at  (OC,  32F)  freezing 

%  rho_US=0 . 0588;  %  For  9000  feet  at  (-6.67C,  20F) 

%  rho_US=0.0601;  %  For  9000  feet  at  (-12.22C,  lOF) 

%  rho_US=0.0614;  %  For  9000  feet  at  (-17.77C,  OF) 

%  rho_US=0 . 0628;  %  For  9000  feet  at  (-23.33C,  -lOF) 

g, 

o 
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%  range  for  rho  US: 

%  from  0.0554  (warm,  high  altitude)  to  0.0731  (cold,  low  altitude) 

g, 

o 

%  rho=l*1.15;  %  atmospheric  density  Kg/m^3  from  Muller  (1997) 

(original) 

%  rho=1.2  at  sea  level  and  20  °C 

g, 

o 

%  Density  of  Air  =  P/ {R*T) 

%  D  =  density,  kg/m3 

%  P  =  pressure,  Pascals  {  multiply  mb  by  100  to  get  Pascals) 

%  R  =  gas  constant  ,  J/ {kg*degK)  =  287.05  for  dry  air 

%  T  =  temperature,  deg  K  =  deg  C  +  273.15 


g, 

o 

rho=rho  US*16.0185;  %  kg/m3  %  conversion 

g, 

o 

mass=70;  %  kg  from  Muller (1997) ...I  want  to  have  other  weights 

area=0.4;  %  half  avg  drag  and  lift  area  {m^2)  from  Muller (1997) 

%  I  determined  this  to  be  1.0  m^2  for  a  6'  tall  skier 
%  I  want  to  vary  this  with  time,  as  a  skier  changes 
%  position  during  flight 
%Choose  a  fall  height  and  takeoff  angle 
del  v=l*2.25;  %  additional  pop  off  of  the  jump 
betaod=25;  %  Jump  angle  in  degrees  %  new  name  for  thod 
%  this  will  vary 

betao=betaod*pi/180; %  takeoff  ramp  angle  in  radians 

sbo=sin (betao) ;  %  sine  of  beta 

cbo=cos (betao) ;  %  cosine  of  beta 

tbo=tan (betao) ;  %  tangent  of  beta 

h=l;  %meters 

%hnom=h;  %  nominal  fall  height  to  taper  to 

%  reference  frame  origin  is  at  takeoff  point 
%  choose  initial  surface  point  at  zenith  of  surface 

g, 

o 

%  The  next  sequence  finds  the  critical  point  for  given  parameters 
a=l; 

b=-2*cbo*sbo; 
c=4  *cbo^2 ; 
d=-2*h*cbo*sbo; 
f=2*h*cbo^2; 

g, 

o 

xc= (c*d-b*f ) / (b^2-a*c)  %  center  (x)  of  ellipse  for  critical  point 
search 

yc= (a*f-b*d) / (b^2-a*c)  %  center  (y)  of  ellipse  for  critical  point 
search 

%  rotation  angle  (alpha) :  from  x-axis  to  major  axis 
%  same  as  from  y-axis  to  minor  axis 
alpha=. 5*atan ( (2*b) / (a-c) ) ;  %  in  radians 
alpha  degrees=alpha*180/pi; 

j=  2*7 (a*f"2) + (c*d"2) - (2*b*d*f) ) ;%numerator 

k= (b^2-a*c) * (-sqrt ( ( (a-c) ^2) + (4*b^2) ) - (a+c) ) ;  %  denominator 

g, 

o 

%  L  =  half  of  semi  minor  axis  length 
L= (sqrt ( j /k) ) ;  %L= (sqrt ( j /k) /2) 

g, 

o 
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critical  pt  x=xc+ {L*sin (alpha) ) 
critical_pt_y=yc- {L*cos (alpha) ) 

g, 

o 

%  Choose  a  set  of  velocities 
for  yL=[5,  2.5,  0,  -2.5,  -5]; 

%for  yL=[5,  0,-5, -10, -15]; 

%yL=-8 . 5 

%for  yL= [ 1 ,  -1, -3 . 5, -7, -8 . 5] ;  %  won't  work  for  xL=40 

%  5  yL ' s  at  x=40,  bring  to  x=10.  Overshoot  with  some  velocity 
%  higher  and  closer  together  that  xL  should  be 
%  xL=40;  look  at  x=10;  all  between  2  and  3. 

%  purpose  is  to  generate  four  safe  slopes,  at  vo=  20,16,12,8 

g, 

o 

%xL=15; 

xL=40;  %xL=50;  %  horizontal  length  of  jump  (m) 

g, 

o 

dx=l;%  dx=.01;  %  spacing  of  slope  calculation  interval 

x=xL : -dx : critical  pt  x+eps;  %  backwards  integration  from  xL  to  crit. 

pt. 

options  =  odeset ( ' MaxStep ' , 1 ) ; 

% [xsf , ysf ] =odel5s (@safesurface, x, yL, options) ; 

[xsf , ysf] =ode45 (@safesurface, x, yL, options)  ; 
xysf= [xsf  ysf ] ; 

xsf=f lipud (xsf ) ;  %  integrate  from  right  to  left 

ysf=f lipud (ysf ) ;  %  backwards  because  integrating 

%  makes  a  difference  in  uses  of  safe  surfaces  later  on,  in  impact 
events 

g, 

o 

xs=xsf;  %  flipped  up  above 
ys=ysf ; 

g, 

o 

slopeshape= [xsf  ysf]  ;  %  calculate  const  fall  ht  surface  %  is  this 

needed? 

dy=diff (ysf) ; 

g, 

o 

dysdxs= [dy ( 1 ) ;  dy] ;  %  note  assumes  dx=l  %  is  this  even  needed? 

figure ( 1 ) ; 

plot (xsf , ysf , ' k ' , ' linewidth ' , 2 . 0 ) ;  %  plot  in  meters 
axis  equal; 
grid  on; 

xlabel (  '  X  (m) ' ) ; 
ylabel (  '  y  (m) ' ) ; 
hold  on; 
end 

g, 

o 

%  plot  trajectories  and  verify  efh 
%  integrate  flight  equations  including  drag  to 
%  check  whether  safe  slope 
%  designed  neglecting  drag  is  still  safe 
heq=l ; 
dx=0 . 01 ; 
tmax=9; 
dt=tmax/50; 

options  =  odeset (' Events ', @impactevents) ; 

for  i=l:l:8;  %  add  pop  by  adding  delta  theta  to  betao 
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vimph=i*5;%  max  speed  40  mph 
vi=vimph*mph2mps ; 

delta  theta=atan (del  v/vi) ;  %  takeoff  angle  increment  after  pop 
thetao=betao  +  delta  theta  ;  %total  initial  jumper  angle  with 
horizontal 

cto=cos (thetao)  ; 
sto=sin (thetao)  ; 
tto=tan (thetao) ; 

vo=(vi^2+del  v^2)^.5;  %calculates  the  jumper  velocity 

%adding  skier  pop 

vTO (i) =vo*mps2mph;  %vo=voV (i) ; 

tspan= [ 0 : dt : tmax] ; 

vox=vo*cto; 

voy=vo*sto; 

ics=[0  vox  0  voy  ]  ; 

[t, y, tE, yE, iE] =ode45 (@ jumpflight, tspan, ics, options) ; 

eventE=[tE  yE  iE] ; 

xtraj=y ( : , 1) ; 

ytraj=y ( ; ,  3)  ; 

tflight (i) =t (end)  ; 

%tflight (i) =tE  (2) ; 

vximp=yE (2 , 2 ) ;  %  impact  is  always  the  second  event  to  occur 

vyimp=yE (2,4)  ; 

lenimp (i) =sqrt (yE (2,1) ^2+yE (2, 3) ^2) ; 

%lenimpf t ( i ) =lenimp/ ft 2m; 

vimp=sqrt (vximp^2+vyimp^2 ) ;  %  impact  velocity 
thimp=atan (  yE (2 , 4 ) /yE (2 , 2 )  ) ;  %  impact  angle 

tpk=tE(l);  %  peak  ht  is  always  first  event 
ypk=interpl (t, ytraj , tpk, ' spline ' ) ; 
xpk=interpl (t, xtraj , tpk, ' spline ' ) ; 
ysurfpk=interpl (xs, ys, xpk,  ' spline ' )  ; 
airht (i) =ypk-ysurfpk; 

%  independently  calculate  slope  angle  at  impact  point  from  xs,  ys 
%dx=0.001;  %meters 
yplusdy=spline (xs, ys, yE (2, 1) +dx) ; 
dydx= (yplusdy-yE (2,3) ) /dx; 

slopeangledeg=atan (dydx) /deg2rad;  %  point  of  this  is  to  calculate 
%  vperp  and  equivalent  fall  ht . 
slopeR ( i ) =slopeangledeg; 

vperp=-vimp*sin (thimp-slopeangledeg*deg2rad)  ; 

heq (i) = (vperp. ^2) / (2*g)  ;  %heq  is  equivalent  fall  ht . 

g, 

o 

figure  (2 )  ; 
hold  on 

subplot (2,1,1); 

plot (lenimp ( i )  ,  heq ( i )  ,  ' ko ' , 0 , 0 ) ; 

xlabel ( ' jump  length  (m) ' ) ; 

ylabel ( '  equivalent  fall  height  (m) ' ) ; 

title(['For  takeoff  angle  =  ' , num2str (betao*180/pi) , ' deg  and 
nominal  fall  height  =',  num2str (hnom) , 'm' ] ) ; 
hold  off 

g, 

o 

hold  on 

subplot (2,1,2); 
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plot { lenimp { i ) , -slopeangledeg,  ' ko ' ) ; 
xlabel { ' jump  length  (m) ' ) ; 
ylabel { '  landing  slope  angle  (deg)  ' )  ; 
hold  off 

figure { 3 ) ; 
hold  on 

subplot (2,1,1) ; 

plot (lenimp ( i ) , tf light ( i )  ,  ' ko ' , 0 , 0 ) ; 
xlabel ( ' jump  length  (m) ' ) ; 
ylabel ( '  flight  time  (s)'); 

title(['For  takeoff  angle  =  ' , num2str (betao*180/pi) , ' deg  and 
nominal  fall  height  =',  num2str (hnom) , 'm' ] ) ; 
hold  off 
hold  on 

subplot (2,1,2) ; 

plot (lenimp ( i ) , airht ( i )  ,  ' ko '  )  ; 
xlabel ( ' jump  length  (m) ' ) ; 
ylabel ( '  jump  air  height  (m) ' ) ; 
hold  off 

figure  ( 1 ) ; 
subplot (1,1,1)  ; 
hold  on 
axis  equal 

plot (xt raj , yt raj , ' k ' ) ; 
xlabel  (  '  X  (m)  ' ) 
ylabel ( '  y  (m) ' ) 

title(['For  takeoff  angle  =  ' , num2str (betao*180/pi) , ' deg  and 
nominal  fall  height  =',  num2str (hnom) , ' Im' ] ) 
hold  off 

end 

takeoff  angle  nominal  fallht= [betao*rad2deg  h] ; 

results= [vTO '  airht'  heq'  lenimp'  slopeR'  tflight']; 

[mm, nnn] =size (results)  ; 

disp ( sprintf ( ' For  takeoff  angle  =  %-5.1f(deg)  and  nominal  fall 
height  =  %-4.1f(m)',  betao*180/pi, hnom) ) 

disp ( sprintf (' vTO (mph)  airht (m)  heq (m)  lenimp (m)  slope (deg) 

tflight ( s )  '  ) ) 

for  1=1 :mm 

disp( sprint f('%-10.4f',  re sults(i,  :)  )  ) 

end 

delxL=xL- lenimp; 
delhimp=h-heq; 


B.  SINGULAR  POINT  LOCATION  CODE 

format  compact; 

rad2deg=180/pi;  %  converts  radians  to  degrees 
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ft2m=0 . 3048;  %  converts  feet  to  meters 

m2ft=l/ft2m;  %  converts  meters  to  feet 

deg2rad=pi/180; %  converts  degrees  to  radians 

mps2mph=2 . 2369362 9;  %  converts  meters/second  to  miles/hour  (changed 
from  2.237) 

mph2mps=l/mps2mph; %  miles  per  hour  to  meters  per  second  conversion 
g=9. 80665;  %  meters  per  second  (changed  from  9.81) 

g, 

o 

%  to  convert  lb/ft3  to  kg/m3,  multiply  by  16.0185 
%  to  convert  kg/m3  to  lb/ft3,  multiply  by  0.0624 

%  From:  http://www.denysschen.com/catalogue/density.asp 

g, 

o 

%  Uncomment  any  rho  US  desired,  depending  on  elevation  of  jump 
%  rho  US  is  measured  in  lb/ft^3  ;  see  conversion 


rho_US=0 .0597; 

g, 

o 


rho=rho  US*16.0185; 

g, 

o 


For  8000  feet  at  (OC, 
%  kg/m3  %  conversion 


32F) 


freezing 


mass=70;  %  kg  from  Muller (1997) ...I  want  to  have  other  weights 

area=0.4;  %  half  avg  drag  and  lift  area  (m^2)  from  Muller (1997) 

%Choose  a  fall  height  and  takeoff  angle 
del  v=l*2.25;  %  additional  pop  off  of  the  jump 
betaod=25;  %  Jump  angle  in  degrees  %  new  name  for  thod 
%  this  will  vary 

betao=betaod*pi/180; %  takeoff  ramp  angle  in  radians 

sbo=sin (betao) ;  %  sine  of  beta 

cbo=cos (betao) ;  %  cosine  of  beta 

tbo=tan (betao) ;  %  tangent  of  beta 

h=l;  %meters 

%hnom=h;  %  nominal  fall  height  to  taper  to 

%  reference  frame  origin  is  at  takeoff  point 
%  choose  initial  surface  point  at  zenith  of  surface 

g, 

o 

%  The  next  sequence  finds  the  critical  point  for  given  parameters 
a=l; 

b=-2*cbo*sbo; 


c=4*cbo^2; 

d=-2*h*cbo*sbo; 

f=2*h*cbo^2; 

g, 

o 

xc= (c*d-b*f ) / (b^2-a*c) ;  %  center  (x)  of  ellipse 
search 

yc= (a*f-b*d) / (b^2-a*c) ;  %  center  (y)  of  ellipse 
search 

%  rotation  angle  (alpha) :  from  x-axis  to  major 
%  same  as  from  y-axis  to  minor  axis 
alpha=. 5*atan ( (2*b) / (a-c) ) ;  %  in  radians 
alpha  degrees=alpha*180/pi; 

j=  2*7 (a*f"2) + (c*d"2) - (2*b*d*f) ) ;%numerator 
k= (b"2-a*c) * (-sqrt ( ( (a-c) "2) + (4*b^2) ) - (a+c) ) ;  % 

g, 

o 

%  L  =  half  of  semi  minor  axis  length 
L= (sqrt ( j /k) ) ;  %L= (sqrt ( j /k) /2) 


for  critical 
for  critical 
axis 


denominator 


point 

point 
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singular  pt  x=xc+ {L*sin (alpha) ) 
singular_pt_y=yc- {L*cos (alpha) ) 


C.  VELOCITY  WITH  DRAG  ALONG  IN-RUN  CODE 

%  rho  US  is  measured  in  lb/ft^3  ;  see  conversion 

rho_US=0 . 05 97 ;  %  For  8000  feet  at  (OC,  32F)  freezing 

rho=rho  US*16.0185;  %  kg/m3  %  conversion 

A=0.5;  %  cross  sectional  area,  m^2 

Cd=0.42;%  Drag  coefficient,  for  a  half  sphere 

g=  9 . 8  0  6  6  5 ; 

lambda_degrees=25 ; 

lambda=lambda  degrees*pi/180; 

muu=0 . 05 ; 

m=40;  %kg 

gg=g*sin (lambda) -muu*g*cos (lambda) ; 
t= [ 0 : 0 . 1 : 2 ] ; 

K=rho*A*Cd; 

L=10;  %  in-run  (meter) 

g, 

o 

a=sqrt (2*m*gg/ (rho*A*Cd) ) ; 

g, 

o 

distance= ( (a^2 ) /gg) *log (cosh (gg/a*t) ) ;  %  distance  for  a  time  t 
velocity=sqrt (2*m*gg/K) *tanh (t*sqrt (K*gg/2*m) ) ;  %  velocity  for  a  time  t 

velocity2=sqrt (2*m*gg/K) *tanh (acosh (exp (L*K/ (2*m) ) ) ) 

%  velocity  for  a  linear  in-run  length  L 


D.  EQUIVALENT  FALL  HEIGHT  CODES 
I.  Tabletop  EFH  Code 

%  EFH  for  tabletop  jump 
%  Andrew  Swedberg 
%  NFS,  April  2010 

%  yt=height  above  tabletop  (m) ,  written  as  negative  value 
%  xt  is  length  of  tabletop  (m) 

%  phi  is  angle  of  landing  slope  (radians),  written  as  negative  value 
%  theta  zero  is  angle  of  takeoff  ramp  (degrees) 

%  EFH  is  equivalent  fall  height 
%  X  is  range  for  length  of  jump,  (m) 

g, 

o 

clear 

format  compact 
g=  9 . 8  0  6  6  5 ; 
x= [ 0 : 1 : 4  0 ] ; 
grid  on 

g, 

o 


yt=0;  %  0,  -1,  -2,  -3;  %  meters  (negative) 
phi_degrees=-30;  %  degrees  (negative) 
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degrees 


theta_zero_degrees=25;  % 

g, 
o 

g, 
o 

phi=phi  degrees* {pi/180) ;  %  converts  to  radians 
theta_zero=theta_zero_degrees* (pi/180) ;  %  converts  to  radians 
cto=cos (theta_zero) ; 
tto=tan (theta_zero)  ; 
sto=sin (theta_zero)  ; 

g, 

o 

h2=-yt; 

hl= ( (x. ^2) *tto"2) . / (4* (x. *tto+h2) ) ; 

EFH=hl+h2; 

g, 

o 

hold  on 
figure  (1) 

plot  (x, EFH, ' k . ' , ' linewidth '  ,2.0)  ; %  axis  equal 

axis ([040015])  ; 

xlabel  {'horizontal  length  of  jump  x  (m) ') 
ylabel  { 'EFH  (m) ' ) ; 

g, 

o 

g, 

o 

g=  9 . 8  0  6  6  5 ; 
x= [ 0 : 1 : 4  0 ] ; 

g, 

o 

xt=l 0 ; 

yt=-l;  %  0,  -1,  -2,  -3;  %  meters  (negative) 
phi_degrees=-30;  %  degrees  (negative) 
theta_zero_degrees=25;  %  degrees 

g, 

o 

phi=phi  degrees* {pi/180) ;  %  converts  to  radians 
theta_zero=theta_zero_degrees* (pi/180) ;  %  converts  to  radians 
cto=cos {theta_zero) ; 
tto=tan {theta_zero)  ; 
sto=sin {theta_zero)  ; 

g, 

o 

h2=-yt; 

hl= { (x. "2) *tto"2) . / (4* (x. *tto+h2) ) ; 

EFH=hl+h2; 

g, 

o 

hold  on 
figure  (1) 

plot  {x,EFH, ' k+ ' , 'linewidth'  ,2.0) 
axis (  [0  40  0  12] )  ; 

xlabel  {'horizontal  length  of  jump  x  (m)  ') 
ylabel  { 'EFH  (m) ' ) ; 

g, 

o 

g, 

o 

g=  9 . 8  0  6  6  5 ; 
x= [ 0 : 1 : 4  0 ] ; 

g, 

o 

xt=10; 

yt=-2;  %  0,  -1,  -2,  -3;  %  meters  (negative) 
phi_degrees=-30;  %  degrees  (negative) 
theta_zero_degrees=25;  %  degrees 

g, 

o 
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phi=phi  degrees* {pi/180) ;  %  converts  to  radians 
theta_zero=theta_zero_degrees* (pi/180) ;  %  converts  to  radians 
cto=cos (theta_zero) ; 
tto=tan (theta_zero)  ; 
sto=sin (theta_zero)  ; 

g, 

o 

h2=-yt; 

hl= ( (x. "2) *tto"2) . / (4* (x. *tto+h2) ) ; 

EFH=hl+h2; 

hold  on 
figure  (1) 

plot  (x,  EFH,  '  k'' '  ,  '  linewidth  '  ,2.0)  ;  %  axis  equal 

axis ( [0  40  0  12] )  ; 

xlabel  ('horizontal  length  of  jump  x  (m) ') 
ylabel  { 'EFH  (m) ' ) ; 


g, 

o 

g=  9 . 8  0  6  6  5 ; 
x= [ 0 : 1 : 4  0 ] ; 

g, 

o 

xt=10; 

yt=-3;  %  0,  -1,  -2,  -3;  %  meters  (negative) 
phi_degrees=-30;  %  degrees  (negative) 
theta_zero_degrees=25;  %  degrees 

g, 

o 

g, 

o 

phi=phi  degrees* {pi/180) ;  %  converts  to  radians 
theta_zero=theta_zero_degrees* {pi/180) ;  %  converts  to  radians 
cto=cos (theta_zero) ; 
tto=tan {theta_zero)  ; 
sto=sin {theta_zero)  ; 

g, 

o 

h2=-yt; 

hl= { (x. "2) *tto"2) . / (4* (x. *tto+h2) ) ; 

EFH=hl+h2; 

hold  on 
figure  (1) 

plot  (x, EFH, ' ko ' , ' linewidth '  ,2.0);  %axis  equal 
axis { [0  40  0  12] ) ; 

xlabel  {'horizontal  length  of  jump  x  (m) ') 
ylabel  { 'EFH  (m) ' ) ; 

legend {'yt=0  m  ','yt=l  m','yt=2  m','yt=3  m'); 

2.  Angled  Linear  Downslope  Landing  Code 

%  EFH  for  tabletop  jump 
%  Andrew  Swedberg 
%  NFS,  April  2010 

%  yt=height  above  tabletop  (m) ,  written  as  negative  value 
%  xt  is  length  of  tabletop  (m) 
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%  phi  is  angle  of  landing  slope  (radians),  written  as  negative  value 
%  theta  zero  is  angle  of  takeoff  ramp  (degrees) 

%  EFH  is  equivalent  fall  height 
%  X  is  range  for  length  of  jump,  (m) 

g, 

o 

clear 

format  compact 
g=  9 . 8  0  6  6  5 ; 
x= [ 0 : 1 : 4  0 ] ; 
grid  on 

g, 

o 

xt=l 0 ; 

yt=0;  %  meters  (negative) 

phi_degrees=-30;  %  degrees  (negative) 

theta_zero_degrees=25;  %  degrees 

phi=phi  degrees* (pi/180) ;  %  converts  to  radians 

theta_zero=theta_zero_degrees* (pi/180) ;  %  converts  to  radians 

cto=cos (theta_zero) ; 

tto=tan (theta_zero)  ; 

sto=sin (theta_zero)  ; 

[rl,  cl,  vl] =find (x<xt) ; 

[r2,  c2,  v2] =find (x>=xt) ; 


yl=vl*yt; 

y2=yt+tan (phi) * (x (find (x>=xt) ) -xt) ; 
y=  [yl  y2]  ; 

g, 

o 

phil=vl*0; 
phi2=v2  *phi ; 
phi=[phil  phi2]; 

vo=sqrt ( (x. ^2) *g. / (2* (x. *tto-y) *cto^2) ) ;  %  init  v  to  get  to  x  &  y 
%  velocity  at  x  and  y 

vl=sqrt (vo . ^2 *cto^2+ (vo . *sto-g*x . /vo . /cto)  .  ^2 )  ;  %original  way  to  find 
velocity 

v2=sqrt ( (vo . ^2 ) -2 *g*y) ;  %  impact  velocity,  same  as  vl 

v3=sqrt ( ( (x. ^2) *g. / (2* (x. *tto-y) *cto^2) ) -2*g*y) ; %  alternate  way;  in  my 
thesis,  in  terms  of  x  and  y 

g, 

o 

theta=atan ( tto-g*x . / (vo . ^2 *cto^2 ) ) ;  %original  way  to  find  theta,  impact 
angle 

theta2=atan (2*y. /x-tto) ;  %  alternate  way--this  matches  the  original 

way 

g, 

o 

V  perp=v3 . *sin (phi-theta) ;  %  perpendicular  velocity 
EFH1=  (v  perp.  ''2)  /  (2*g)  ; 

figure  (1) 
hold  on 

plot  (x, EFHl , ' k . ' , ' linewidth '  ,2.0) 
xlabel  ('horizontal  length  of  jump  x  (m) ') 
ylabel  ( 'EFH  (m) ' ) ; 
axis  (  [0  4  0  0  12] ); 
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g=  9 . 8  0  6  6  5 ; 
x= [ 0 : 1 : 4  0 ] ; 

g, 

o 

xt=10; 

yt=-l;  %  meters  (negative) 

phi_degrees=-30;  %  degrees  (negative) 

theta_zero_degrees=25;  %  degrees 

phi=phi  degrees* (pi/180) ;  %  converts  to  radians 

theta_zero=theta_zero_degrees* (pi/180) ;  %  converts  to  radians 

cto=cos (theta_zero) ; 

tto=tan (theta_zero)  ; 

sto=sin (theta_zero)  ; 

[rl, cl, vl] =find (x<xt) ; 

[ r2 , c2 , v2 ] =f ind (x>=xt ) ; 


yl=vl*yt; 

y2=yt+tan (phi) * (x (find (x>=xt) ) -xt) ; 
y=  [yl  y2]  ; 

g, 

o 

phil=vl*0; 
phi2=v2  *phi ; 
phi=[phil  phi2]; 

vo=sqrt ( (x. ^2) *g. / (2* (x. *tto-y) *cto^2) ) ;  %  init  v  to  get  to  x  &  y 
%  velocity  at  x  and  y 

vl=sqrt (vo . ^2 *cto^2+ (vo . *sto-g*x . /vo . /cto) . ^2 ) ;  %original  way  to  find 
velocity 

v2=sqrt ( (vo . ^2 ) -2 *g*y) ;  %  impact  velocity,  same  as  vl 

v3=sqrt ( ( (x. ^2) *g. / (2* (x. *tto-y) *cto^2) ) -2*g*y) ; %  alternate  way;  in  my 
thesis,  in  terms  of  x  and  y 

g, 

o 

theta=atan ( tto-g*x . / (vo . ^2 *cto^2 ) ) ;  %original  way  to  find  theta,  impact 
angle 

theta2=atan (2*y. /x-tto) ;  %  alternate  way--this  matches  the  original 

way 

g, 

o 

V  perp=v3 . *sin (phi-theta) ;  %  perpendicular  velocity 
EFH1=  (v  perp.  ''2)  /  (2*g)  ; 

figure  (1) 

plot  (x, EFHl , ' k+ ' , ' linewidth '  ,2.0) 
hold  on 

xlabel  ('horizontal  length  of  jump  x  (m) ') 
ylabel  ( 'EFH  (m) ' ) ; 
axis (  [0  40  0  12] )  ; 


oooooooooooooooooooooooooooooooooooooo 

g=9 . 80665; 
x= [ 0 : 1 : 4  0 ] ; 

g, 

o 
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xt=10; 

yt=-2;  %  meters  (negative) 

phi_degrees=-30;  %  degrees  (negative) 

theta_zero_degrees=25;  %  degrees 

phi=phi  degrees* (pi/180) ;  %  converts  to  radians 

theta_zero=theta_zero_degrees* (pi/180) ;  %  converts  to  radians 

cto=cos (theta_zero) ; 

tto=tan (theta_zero)  ; 

sto=sin (theta_zero)  ; 

[ rl , cl , vl ] =f ind (x<xt )  ; 

[ r2  ,  c2 , v2 ] =f ind (x>=xt ) ; 


yl=vl*yt; 

y2=yt+tan (phi) * (x (find (x>=xt) ) -xt) ; 
y=  [yl  y2]  ; 

g, 

o 

phil=vl*0; 
phi2=v2  *phi ; 
phi=[phil  phi2]; 

vo=sqrt ( (x. ^2) *g. / (2* (x. *tto-y) *cto^2) ) ;  %  init  v  to  get  to  x  &  y 
%  velocity  at  x  and  y 

vl=sqrt (vo . ^2 *cto^2+ (vo . *sto-g*x . /vo . /cto)  .  ^2 )  ;  %original  way  to  find 
velocity 

v2=sqrt ( (vo . ^2 ) -2 *g*y) ;  %  impact  velocity,  same  as  vl 

v3=sqrt ( ( (x. ^2) *g. / (2* (x. *tto-y) *cto^2) ) -2*g*y) ; %  alternate  way;  in  my 
thesis,  in  terms  of  x  and  y 

g, 

o 

theta=atan ( tto-g*x . / (vo . ^2 *cto^2 ) ) ;  %original  way  to  find  theta,  impact 
angle 

theta2=atan (2*y. /x-tto) ;  %  alternate  way--this  matches  the  original 

way 

g, 

o 

V  perp=v3 . *sin (phi-theta) ;  %  perpendicular  velocity 
EFH1=  (v  perp.  ''2)  /  (2*g)  ; 

figure  (1) 

plot  (x, EFHl , ' ' linewidth '  ,2.0) 
hold  on 

xlabel  ('horizontal  length  of  jump  x  (m) ') 
ylabel  ( 'EFH  (m) ' ) ; 
axis (  [0  4  0  0  12] )  ; 


9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9' 

oooooooooooooooooooooooooooooooooooooo 

g=  9 . 8  0  6  6  5 ; 
x= [ 0 : 1 : 4  0 ] ; 

g, 

o 

xt=l 0 ; 

yt=-3;  %  meters  (negative) 

phi_degrees=-30;  %  degrees  (negative) 

theta_zero_degrees=25;  %  degrees 

phi=phi  degrees* (pi/180) ;  %  converts  to  radians 

theta_zero=theta_zero_degrees* (pi/180) ;  %  converts  to  radians 
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cto=cos {theta_zero)  ; 
tto=tan {theta_zero)  ; 
sto=sin {theta_zero)  ; 

[rl, cl, vl] =find {x<xt) ; 

[ r2 , c2 , v2 ] =f ind {x>=xt ) ; 


yl=vl*yt; 

y2=yt+tan (phi) * (x (find {x>=xt) ) -xt) ; 
y=  [yl  y2]  ; 

g, 

o 

phil=vl*0; 
phi2=v2  *phi ; 
phi=[phil  phi2]; 

vo=sqrt { (x. ^2) *g. / (2* (x. *tto-y) *cto^2) ) ;  %  init  v  to  get  to  x  &  y 
%  velocity  at  x  and  y 

vl=sqrt (vo . ^2 *cto^2+ (vo . *sto-g*x . /vo . /cto)  .  ^2 )  ;  %original  way  to  find 
velocity 

v2=sqrt { (vo . ^2 ) -2 *g*y) ;  %  impact  velocity,  same  as  vl 

v3=sqrt { { (x. ^2) *g. / (2* (x. *tto-y) *cto^2) ) -2*g*y) ; %  alternate  way;  in  my 
thesis,  in  terms  of  x  and  y 

g, 

o 

theta=atan { tto-g*x . / (vo . ^2 *cto^2 ) ) ;  %original  way  to  find  theta,  impact 
angle 

theta2=atan {2*y. /x-tto) ;  %  alternate  way--this  matches  the  original 

way 

g, 

o 

V  perp=v3 . *sin (phi-theta) ;  %  perpendicular  velocity 
EFH1= (v  perp. ^2) / {2*g) ; 

figure  (1) 
hold  on 

plot  (x, EFHl , ' ko ' , ' linewidth '  ,2.0) 
xlabel  {'horizontal  length  of  jump  x  (m) ') 
ylabel  { 'EFH  (m) ' ) ; 
axis { [0  40  0  12] ) ; 

legend {'yt=0  m  ','yt=l  m','yt=2  m','yt=3  m'); 


E.  CIRCULAR  TRANSITION  CODES 

I.  Circular  Transition  Before  Takeoff 

function  [T,Y]=  before  takeoff 

format  compact 

%clf 

global  g  muu  lambda  beta  gamma  Z  radius  min 
lambda_degrees=35 ; 
beta_degrees=10 ; 

lambda=lambda_degrees* {pi/180) ;  %  radians 
beta=beta_degrees* {pi/180) ;  %  radians 
muu=0.05;  %  coefficient  of  friction  of  skis  on  snow 
gamma=1.5;  %  This  means  1.5  times  the  force  of  gravity 
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g=9. 80665;  %  gravity 

%  Compute  the  initial  velocity  at  the  beginning  of  the  transition 
Z=10;  %  10  meter  in-run  part 

%  Inputs: 

%  radius  min  =  minimum  radius  of  circle  in  meters 
%  lambda  =  in  degrees 
%  Returns : 

%  T  =  Time  vector 

%  Y  =  Vector  of  Velocities  and  Epsilons 

vti=sqrt {2*g*Z* (sin (lambda) -muu*cos (lambda) ) ) 

%vti=sqrt (2*g. *Z . * (cos (lambda) -muu . *sin (lambda) ) ) 
radius  min= (vti . ^2 ) / (gamma*g) 

%  pause 

%  vo=velocity  at  the  end  of  the  in-run  =  beginning  of  the  transition 

%  Set  up  the  stopping  criteria 
options=odeset ( ' Events ' ,  @circular_events ) ; 

%  Set  up  the  initial  conditions 
tspan= [0 : .01:5] ; 

IC=[vti;-l ambda ]  ; 

%initial  conditions  for  the  transition  part 

%  velocity  it  at  vo,  angle  (epsilon)  begins  at  negative  lambda 
%  Set  up  the  time  window  to  solve  the  ODE. . .  may  not  take  this  long. . . 
%  better  longer  than  shorter 

%  Make  the  ode45  solver  call  --  note  that  all  parameters  are  passed. 

%  Stopping  criteria  is  when  epsilon  reaches  +beta 
[T, Y] =ode4  5 ( @partl , tspan,  IC,  options )  ; 

Addl  Accel= (Y (:, 1 ). ^2 ) /radius  min; 

g, 

o 

time  at  beta=T(end);  %  Final  time  where  epsilon  reaches  +beta 
%beta-Y (end, 2 ) ;  %  Value  of  epsilon  -  beta  (should  be  zero) 

g, 

o 

figure ( 1 ) 

plot (T, Y ( : , 1 ) , ' k ' , ' linewidth ' ,  2.0)  %  Velocity 
%title (' velocity  vs.  time  in  a  circular  transition') 
xlabel ( ' time,  (s) ') 
ylabel (' velocity,  (m/s)') 

g, 

o 

figure (2 ) 

plot(T, (g*cos(Y(:,2)) +Addl_Accel ) , ' k ' ,  ' linewidth ',2.0) 

%title (' Normal  acceleration  vs.  time  in  a  circular  transition') 
xlabel (' time,  (s) ') 
ylabel (' Normal  Force,  (m/s^2) ') 

g, 
o 

figure ( 3 ) 

Y  (  :  ,  2 )  ; 
hold  on 

for  i=[-2:0.01:0] ; 

plot  ( i , cos ( lambda) ,' k ' ,  ' linewidth ', 1 . 0 ) 

end 

g, 
o 

partial_plot_vector= [T-1 .9568] ; 
plot  (T,  ( (g*cos (Y ( : , 2) ) +Addl_Accel) /g) ,  ' k.  ' ,  ' 

no 


linewidth ',1.0) 


for  j= [T (end) : 0 . 01 : 3] ; 

plot  { j , cos { lambda) ,  'k',  'linewidth',  1.0) 

end 

g, 

o 

axis  { [-2  3  0  3] ) 

%title {' Ratio  of  normal  accleration  to  g  vs.  time  in  a  circular 
transition ' ) 
xlabel { ' time,  (s) ') 

ylabel {' Ratio  of  Normal  force  to  g') 

g, 

o 

function  dvdtl=partl (t, v, muu, g, radius  min, lambda, beta) 
global  g  muu  lambda  beta  radius  min 

vv=v { 1 ) ; 
epsilon=v (2 ) ; 

V  prime=-g*sin (epsilon) -muu* { { {vv^2 ) / radius  min) +g*cos (epsilon) ) ; 
epsilon  prime=vv/radius  min; 
dvdtl=[v  prime;  epsilon  prime]; 

%  Output  (Y)  is  a  matrix  with  two  columns 

%  Column  1  is  the  velocity  during  the  transition  part 

%  Column  2  is  the  range  of  epsilon,  that  goes  from  -lambda  to  +beta 

%  Need  to  use  the  last  value  of  the  velocity  for  the  start  of  the 

%  integration  to  get  the  value  of  vo . 

2.  Circular  Transition  Events 

function [value,  is terminal, direction] =events (T, Y) 
global  lambda  beta 
value ( 1 ) =Y (2 ) -beta; 
isterminal ( 1 ) =1 ; 
direction (1) =0; 
return 

F.  CLOTHOIDAL  TRANSITION  CODES 
I.  Radius  Required  for  Given  y 


clear 

velocity= [5:1:15]  ; 
gamma=l . 5 ; 

g=  9 . 8  0  6  6  5 ; 

radius= (velocity. ^2) / (gamma*g) ; 

plot  (velocity, radius ,' k+ linewidth ', 2 . 0 ) ;  axis  equal 
title ('plot  of  radius  of  transition  with  skier  velocity') 
xlabel (' velocity  (m/s)  of  skier  at  beginning  of  transition') 
ylabel (' radius  required  (m)  for  the  given  value  of  gamma') 
hold  on 

velocity= [5:1:15] ; 

gamma=l .75;  %gamma= [1.2:0. 2:2.0] 

g=  9 . 8  0  6  6  5 ; 
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radius= (velocity. ^2) / {gamma*g) ; 

plot  (velocity,  radius k'' linewidth 2 . 0 )  ;  axis  equal 
hold  on 

velocity= [5:1:15] ; 

gamma=2 . 0 ;  %gamma= [1.2:0. 2:2.0] 

g=  9 . 8  0  6  6  5 ; 

radius= (velocity. ^2) / (gamma*g) ; 

plot  (velocity, radius ,  ' ko ' , ' linewidth ' , 2 . 0 ) ;  axis  equal 

legend  (  '  gainma=l .  5  '  ,  '  gainma=l  .75',  '  gainma=2 . 0  '  )  ;  % 

, ' gamma=l . 8 ' , ' gamma=2 . 0 ' ) ; 

%  %  As  the  value  of  gamma  varies  from  1.2  to  2.0  (easier  to  harder  on 
the 

%  %  human  body),  the  radius  required  decreases. 

g, 

o 

%  %  For  lower  velocities,  the  value  of  gamma  does  not  affect  the  radius 
%  %  required  as  much  as  at  high  velocities. 


2.  Clothoid  Transition  Before  Takeoff 

function  [T,Y]  =  clothoid  before  takeoff 
global  g  A  squared  muu  lambda  beta  gamma  Z  LI  L2 

g, 

o 

format  compact 
%clf 

lambda_degrees=35 ; 
beta_degrees=10 ; 

Z=10;  %  in  meters:  in-run  portion 
muu=0.05;  %  coefficient  of  friction 

gamma=1.5;  %  This  means  1.5  times  the  force  of  gravity 
g=9. 80665;  %  gravity 

g, 

o 

lambda=lambda_degrees* (pi/180) ;  %  radians 
beta=beta_degrees* (pi/180) ;  %  radians 
vti=sqrt (2*g*Z* (sin (lambda) -muu*cos (lambda) ) ) ; 

%vti2=sqrt (2*g*Z* (cos (lambda) -muu*sin (lambda) ) ) 

g, 

o 

%  This  is  the  velocity  at  the  beginning  of  the  transition:  where  the 
%  linear  slope  (infinite  radius)  meets  the  clothoidal  shape 

g, 

o 

radius  min= (vti . ^2 ) / (gamma*g) ; 

g, 

o 

%  The  radius  min  is  the  shortest  that  the  "circular"  transition 
%  will  be.  Every  other  radius  length  will  be  longer  than  that,  as  this 
%  will  ensure  the  g-force  felt  by  the  skier  is  always  less  than  a 
%  desired  gamma  value. 

g, 

o 

A  squared= ( radius  min^2 ) * ( lambda+beta) ; 

A=sqrt (A  squared) ;  %  A  is  the  flatness  of  the  curve 
L1=A  squared/radius  min; 

L2=A  squared/radius  min; 
clothoid  length=A*sqrt (lambda+beta) ; 
total  transition  length=2 *clothoid  length 
constantl= (lambda+beta) / (2*L1^2) ; 
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constant2= { lambda+beta) / {2*L2^2)  ; 

g, 

o 

options=odeset {' Events @clothoid_events2 ) ;  %  looks  for  this  event 
tspanl= [ 0 : .05:4] ; 

icl= [vti  0];  %  initial  conditions  for  the  first  clothoid 

[T, Y, TE, YE, IE] =ode45 { @clothoid_odel , tspanl , id , options) ; 

TY= [T  Y] ; 

rl=A_squared . / { Y {  :  ,  2  )  )  ; 
size (rl) ; 

g, 

o 

%  this  is  the  radius--goes  from  infinity  to  radius  min 
%Addl_accel_l= { (Y { : , 2 ) . * (Y { : , 1 ) . ^2 ) ) /A_squared) 

Centrepital_accel_l= { Y { : , 1 ) . ^2 ) . /rl ; 

g, 

o 

%  time  span  is  from  0  to  4  seconds;  it  really  stops  before  this  point 
%  initial  velocity  condition:  vti  (velocity  at  transition) 

%  initial  length  of  clothoid  condition:  0 

g, 

o 

figure { 1 ) 

plot{T,  Y { : , 1 ) , ' k ' , ' linewidth ' ,  2.0); 

%title {' velocity  vs.  time  in  a  clothoid  transition') 
ylabel {' velocity,  (m/s)') 
xlabel { ' time,  (s) ') 

%  The  stopping  criteria  will  be  when  the  clothoid  length  =  s. 

%  For  the  second  clothoid,  the  initial  conditions  will  be 
%  the  velocity  at  the  end  of  the  first  clothoid 

g, 

o 

options=odeset {' Events ', @second_clothoid_events2 ) ;  %  look  for  this 
event 

%  [TT,  YY,  TTE,  YYE,  HE]  = 

ode 4 5 { @clothoid_ode2 , [T (end) ; 6] , [ Y (end) ; L2 ] , options ) 
tspan2= [T (end) : . 05 : 6] ; 

ic2= [Y (end, 1) , Li] ;  %  initial  conditions  for  the  second  clothoid 

g, 

o 

[TT, YY, TTE, YYE, HE]  =  ode45 { @clothoid_ode2 , tspan2, ic2, options) ; 

TTYY=[TT  YY] ; 

g, 

o 

r2=A_squared . / { YY { :  ,  2  )  )  ; 
size {r2) ; 

Centrepital_accel_2= { YY { : , 1 ) . ^2 ) . /r2 ; 

%Addl_accel_22= { (YY { : , 2 ) . * (YY { : , 1 ) . ^2 ) ) /A_squared) 

%Addl  accel  22-Addl  accel  2 
%pause 

g, 

o 

%  time  span  is  from  T (end)  to  6  seconds  (should  stop  before  this  point) 
%  initial  velocity  condition:  velocity  at  point  of  clothoid 
intersection 

%  initial  length  of  clothoid  condition:  Li 

g, 

o 

hold  on 

plot(TT,  YY (:, 1 ),' k linewidth '  ,  2.0); 
ylabel (' velocity,  (m/s)') 
xlabel (' time,  (s) ') 

g, 

o 
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sl=Y{;,2);  %  column  of  s,  the  position  along  the  first  clothoid 
s2=YY{:,2);  %  column  of  s,  the  position  along  the  second  clothoid 
%  starts  at  end 
size  (si)  ; 
size  {s2)  ; 

%pause 

zetal=constantl*  (si  .  ^2)  ;  %  turning  angle  along  first  clothoid 
zeta2=constantl *  (s2  .''2)  ; 

%zeta22=f lipud {constant2*  {s2 . ^2) )  ; 

%pause 

%zeta2=zetal (end) +zeta22 
%pause 

thetal=-lambda+zetal ;  %  angle  of  path  from  horizontal  on  first  clothoid 

wuns=ones (size {s2) ) ; 

zeta2  =  zeta2-wuns*zetal (end)  ; 

%theta2= (thetal (end) *wuns  +  zeta2)  ; 

theta2= (thetal (end) *wuns-zeta2 ) ; %  angle  of  path  from  horizontal: 
second  clothoid 
zeta2=-zeta2 ; 
zetal ; 

zeta_degreesl=zetal* 180/pi 
zeta_degrees2=zeta2* 180/pi 
%  thetal 
%  zeta2 
%  theta2 

g, 

o 

figure (2 ) 
hold  on 

g, 

o 

normal  force  l=g*cos (thetal ) iCentrepital  accel  1;  %  normalized  force 
(mass  specific) ,  divided  by  m. 

normal  force  2=g*cos {theta2 ) iCentrepital  accel  2;  %  second  term  is 
normal  acceleration 

plot (T, normal  force  1,  ' k ' , ' linewidth ' ,  2.0)  ; 

plot (TT, normal  force  2 k linewidth ' ,  2.0)  ; 

%title {' Normalized  mass  specific  force  vs.  time  in  a  clothoid 
transition ' ) 

ylabel {' Force,  {m/s^2)') 
xlabel { ' time,  (s) ') 

g, 

o 

%f igure { 3 ) 
hold  on 

for  i=[-2:0.01:0] ; 

plot  (i, cos (lambda) ,' k' ,  ' linewidth ', 1 . 0 ) 

end 

for  j= [TT (end) : 0 . 01 : 3] ; 

plot  ( j , cos ( lambda) ,  'k',  'linewidth',  1.0) 

end 

g, 

o 

hold  on 

axis  ( [-2  3  0  3] ) 

normal  force  in  g  1= (g*cos (thetal ) +Centrepital  accel  l)/g;  %normal 
force  in  g ' s 

normal  force  in  g  2= (g*cos (theta2 ) iCentrepital  accel  2)/g; 
plot (T, normal  force  in  g  1,'k',  ' linewidth ', 1 . 0 ) 
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linewidth ',1.0) 


plot (TT, normal  force  in  g  2 , ' k ' , 

g, 
o 

%title {' Normal  force  (g)  vs.  time  in  a  clothoid  transition') 
xlabel { ' time,  (s) ') 

ylabel {' Ratio  of  Normalized  mass  specific  force  to  g') 

function  dvdtl  =  clothoid  odel (t, v, lambda, A  squared, g, muu) ; 
global  g  A  squared  muu  lambda  beta  gamma  Z  LI  constant 
vv=v { 1 ) ; 
s=v  (2)  ; 

V  prime=-g*sin { (s . ^2) / {2*A  squared) -lambda) - 

muu*  {  (s  .  /A_squared)  *  (vv.  ''2)  +g*cos  {  (s  .  ^2)  /  {2*A_squared)  -lambda)  )  ; 
s  prime=vv; 

dvdtl=[v  prime; s  prime]; 

g, 

o 

%  the  first  column  of  the  dvdtl  vector  is  v  prime  (s) .  The  second 
column  is  v. 

g, 

o 

function  dvdt2  =  clothoid_ode2 (t, v, beta, A_squared, g, muu) ; 
global  g  A  squared  muu  beta  lambda  gamma  Z  LI 
vv=v { 1 ) ; 
s=v  (2)  ; 

v_prime=g*sin { (s . ^2) / {2*A_squared) -beta) - 

muu*  {  (s  .  /A_squared)  *  (vv.  ''2)  +g*cos  {  (s  .  ^2)  /  {2*A_squared)  -beta)  )  ; 

%  originally  had  a  "-g",  and  positive  betas 
%  v_prime=-g*sin { (s . ^2) / {2*A_squared) +beta) - 

muu*  {  (s  .  /A_squared)  *  (vv.  ''2)  +g*cos  {  (s  .  ^2)  /  {2*A_squared)  +beta)  )  ; 
s  prime=-vv; 

dvdt2= [v  prime; s  prime]; 


3.  Clothoid  Plot 

dt  =  .01  ; 

t= {-2*pi) :dt: {2*pi)  ; 

X  =  cumtrapz { sin { t . ^2 ) )  *  dt  ; 

y  =  cumtrapz (cos (t . ^2 ) )  *  dt  ; 

plot (x, y, ' k ' ,  'linewidth',  2.0)  ;  axis  equal 

axis { [-1  2  -1  2] ) 

4.  Complete  Clothoidal  Transition  Side  View 

%  This  program  will  find  the  shape  of  a  clothoid  transition 

%  for  the  in-run  of  a  ski  jump 

clear 

%  Enter  a  lambda  (existing  slope)  degree 
%  and  a  beta  (desired  take-off  ramp/jump  degree) . 

lambda_degrees=2  0 ; 
beta_degrees=2  0 ; 

g, 

o 

%  si  =  left  side  clothoid  length  at  any  point  (downhill  part) 
%  s2  =  right  side  clothoid  at  any  point  length  (uphill  part) 

%  LI  =  left  side  (longest)  clothoid  length  (downhill  part) 
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%  L2  =  right  side  (longest)  clothoid  length  (uphill  part) 

g, 

o 

%  enter  a  value  for  the  length  of  slope  before  entering  the  transition 

g, 

o 

Z=10;  %  in  meters:  in-run  portion 
mu=0.05;  %  coefficient  of  friction 

gamma=1.5;  %  This  means  1.5  times  the  force  of  gravity 
g=9. 80665;  %  gravity 

g, 

o 

lambda=lambda_degrees* (pi/180) ;  %  radians 
beta=beta_degrees* (pi/180) ;  %  radians 
%  total_turn= ( lambda+beta)  ; 

%  lambda2=total  turn/2; 

%  beta2=total_turn/2 ; 
rotation= (lambda-beta) /2 ; 

g, 

o 

vti=sqrt (2*g. *Z . * (cos (lambda) -mu. *sin (lambda) ) ) ; 

%  This  is  the  velocity  at  the  beginning  of  the  transition:  where  the 
%  linear  slope  (infinite  radius)  meets  the  clothoidal  shape 

g, 

o 

radius  min= (vti . ^2 ) / (gamma*g) 

g, 

o 

%  The  radius  min  is  the  shortest  that  the  "circular"  transition 
%  will  be.  Every  other  radius  length  will  be  longer  than  that,  as  this 
%  will  ensure  the  g-force  felt  by  the  skier  is  always  less  than  a 
%  desired  gamma  value. 

g, 

o 

A  squared= ( radius  min^2 )*(( lambda+beta) ) ;  %  divide  this  by  2??? 

A=sqrt (A_squared) ; 

clothoid  length=A*sqrt (lambda+beta) 
total  transition  length=2 *clothoid  length 

g, 

o 

s= [clothoid_length : -0.1:0]; 

%R= [radius_min: 0.5:1000] ; 

g, 

o 

%  varies  the  value  of  R  from  radius  min  to  a  high  value 
%  (infinite  slope  for  a  linear  portion) 

g, 

o 

%  clothoid:  curve  such  that  the  reciprocal  of  the  radius  varies 
linearly 

g, 

o 

%s=(A^2) ./R;  %  relationship  for  a  clothoid;  used  for  the  range  of  R 

g, 

o 

theta= (max (s) ) / (2*radius_min) ; 

%  since  theta=L/ (2*R) 

g, 

o 

Ll=2*lambda*radius  min;  %  length  of  left  clothoid 
L2=2*beta*radius  min;  %  length  of  right  clothoid 

g, 

o 

Xl=s-(s."5) ./(40.*A."4)+(s."9) . / (3456*A. "8) ; 

%  X  coordinate  for  the  clothoid  shape 

Yl= (s. "3) / (6*A. "2) - (s. "7) / (336*A. "6) + (s . "11) / (42240*A. "10) ; 

%  y  coordinate  for  the  clothoid  shape 

g, 

o 

figure  ( 1 )  ; 
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hold  on 

%plot  of  original  clothoid  (not  rotated) 

%plot (Xl , Y1 , ' r ' , ' linewidth ' , 3 . 0 ) ;  axis  equal;  %X1  is  black 

g, 

o 

%plot  (0, radius  min,  ' x ' ) ;  %  this  'x'  marks  the  center  of  the  circle 


g, 

o 


X2=X1-X1 (1) ; 

%X2=X1+X1 (end) -XI  { 1 )  ; 

Y2=Y1-Y1  (1); 

%plot {X2 , Y2 b linewidth 3 . 0 ) ;  axis  equal;  %X2  is  blue 
%  this  shifts  the  original  clothoid  so  that  the  origin  is  correct 

g, 

o 

X3= (cos (theta) *X2) + (sin (theta) *Y2) ; 

Y3= (-sin (theta) *X2) + (cos (theta) *Y2) ; 

%plot (X3 , Y3 r linewidth 2 . 0 ) ;  axis  equal;  %X3  is  red 
%  this  rotates  the  clothoid  to  theta 

g, 

o 


X4=X3; 

Y4=Y3-radius  min; 

%plot (X4 , Y4 r linewidth 2 . 0 ) ;  axis  equal;  %X4  is  thin  red 
%  this  shifts  the  y-axis  so  that  it  can  be  rotated  about  the  new  origin 

g, 

o 


X5=-X4+2* (X4(l)); 

Y5=Y4; 

%plot (X5 , Y5 g linewidth 2 . 0 ) ;  axis  equal;  %X5 
%this  does  the  same  for  the  right  side 

g, 

o 

X6= (cos (rotation) *X4) + (sin (rotation) *Y4) ; 

Y6= (-sin (rotation) *X4)  +  (cos (rotation) *Y4)  ; 

%plot (X6, Y6, ' c linewidth 2 . 0) ;  axis  equal;  %X6 
%  this  rotates  the  left  side  at  the  correct  angle 

g, 

o 

X7= (cos (rotation) *X5) + (sin (rotation) *Y5) ; 

Y7= (-sin (rotation) *X5) + (cos (rotation) *Y5) ; 

%plot (X7 , Y7 y linewidth 2 . 0 ) ;  axis  equal;  %X7 
%  this  rotates  the  left  side  at  the  correct  angle 

g, 

o 


is  thin  green 


is  cyan 
( lambda2 ) 


is  yellow 
( lambda2 ) 


X8=X6-X6  (1) ; 

Y8=Y6-Y6  (1)  ; 

%plot (X8 , Y8 , colors ( ii ),' linewidth 2 . 0 ) ;  axis  equal;  %X8  is  black 
plot (X8 , Y8 k linewidth 2 . 0 ) ;  axis  equal;  %X8  is  black 
%  this  is  the  plot  of  the  left  clothoid  at  the  correct  angle 

g, 

o 


X9=X7-X7 (1) ; 

Y9=Y7-Y7  (1); 

%plot (X9, Y9, colors (ii) linewidth 2 . 0) ;  axis  equal;%X9  is  black 
plot (X9, Y9, ' k linewidth 2 . 0 ) ;  axis  equal;%X9  is  black 
%title ( ' Plot  of  X  and  y  coordinates  of  clothoid') 
xlabel ( ' X  (m) ' ) 
ylabel ( ' y  (m) ' ) 

%legend ( ' lambda=10 ' , ' lambda=20 ' , ' lambda=30 ' ) ; 

%  this  is  the  plot  of  the  right  clothoid  at  the  correct  angle 
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G.  IN-RUN  WITH  AIR  DRAG  CODE 

m=2  0 ;  %  kg 

%  rho  US  is  measured  in  lb/ft^3  ;  see  conversion 

rho_US=0 . 05 97 ;  %  For  8000  feet  at  (OC,  32F)  freezing 

rho=rho  US*16.0185;  %  kg/m3  %  conversion 

A=0.5;  %  cross  sectional  area,  m^2 

Cd=0.42;%  Drag  coefficient,  for  a  half  sphere 

g=  9 . 8  0  6  6  5 ; 

lambda_degrees=25 ; 

lambda=lambda  degrees*pi/180; 

muu=0 . 05 ; 

gg=g*sin (lambda) -muu*g*cos (lambda) ; 
t= [ 0 : 0 . 1 : 2 ] ; 

K=rho*A*Cd; 

g, 

o 

a=sqrt (2*m*gg/ (rho*A*Cd) ) ; 

g, 

o 

distance= ( (a^2) /gg) *log(cosh(gg/a*t) ) 
velocity=sqrt (2*m*gg/K) *tanh (t*sqrt (K*gg/2*m) ) 
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